AD- A 194  521 


FS-87-193 
SEPTEMBER  19S7 


IIILtlii, 


INFRARED  SHIP  DETECTION  AT  LOW 
SIGNAL-TONOISE 


DTIC 


ELECTE 
APR  2  5  1988 

H 


D 


FLEET  SYSTEMS  DEPARTMENT 

THE  JOHNS  HOPKINS  UNIVERSITY  ■  APPLIED  PHYSICS  LABORATORY 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  ith§cm$ 
Distribution  Unllmltsd 


88 


4  22 


021  t 

i 

RmwsTOnununBnfli 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

lb.  RESTRICTIVE  MARKINGS  1 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 

2b.  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

FS-87-193 

S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 

6b  OFFICE  SYMBOL 
(If  applicable) 

7a.  NAME  OF  MONITORING  ORGANIZATION 

6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Johns  Hopkins  Road 

Laurel,  MD  20707 

7b.  ADDRESS  {City,  State,  and  ZIP  Code) 

8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

8b.  OFFICE  SYMBOL 

Of  applicable) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

8c  ADDRESS  (City,  State,  and  ZIP  Cod*) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 
ELEMENT  NO.  I  NO. 


1 1 .  TITLE  ( Include  Security  Classification) 

Infrared  Ship  Detection  at  Low  Slgnal-to-Noise 


12.  PERSONAL  AUTHOR(S) 

R.  A.  Steinberg  and  J.  J.  fctvera  _ 


13a.  TYPE  OF  REPORT  Il3b.  TIME  COVERED 

Interim  I  from _ TO 


16.  SUPPLEMENTARY  NOTATION 


14.  DATE  OF  REPORT  (Year.  Month,  Day)  |15.  PAGE  COUNT 

1987  September  I  81 


COSATI  CODES 


GROUP  SUBGROUP 


18.  SUBJECT  TERMS  ( Continue  on  reverse  If  necessary  and  Identify  by  block  number) 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  multi -resolution  signal  processing  algorithm  ("MRSI")  was  devised  to 
maximize  the  detection  range  of  imaging  infrared  seekers  against  ship  targets.  Image- 
based  simulations  performed  on  an  IBM  PC  AT  verify  that  the  new  algorithm  provides  a 
signal-to-noise  enhancement  (relative  to  hot-spot  detection)  given  by  ,  where  Npjj£ 

is  the  number  of  pixels  occupied  by  the  target  in  the  image ,  and 

»j-0.81.  The  processor  efficiency  is  nearly  invariant  with  respect  to  sensor/ship  range 
(i.e. , scale) ,  viewing  aspect,  and  ship  class.  Probabilities  of  detection,  determined  via 
2500  image-based  Monte  Carlo  simulations,  agree  with  well-known  analogous  results  of 
psychophysical  experiments  performed  with  human  subjects.  Regardless  of  whether  final 
target  classification  is  performed  on-board  the  missile  via  a  computer  algorithm,  or  off- 
board  by  a  person  inspecting  imagery  telemetered  from  the  missile,  MRSI  substantially 
improves  predicted  system  performance. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

O UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  DTIC  USERS  _ _ 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b.  TELEPHONE  (Indude  Area  Code)  22c.  OFFICE  SYMBOL 

Richard  Steinberg  1(301)  953-6527 _ Croup  El 


OD  FORM  1473, 84  MAR  *3  APR  edition  may  be  uied  until  exhausted.  SECURITY  CLASSIFICATION  OF  THIS  PAGE 


All  other  editions  are  obsolete. 


FS- 87*193 
SEPTEMBER  1987 


INFRARED  SHIP  DETECTION  AT  LOW 
SIGNAL-TO-NOISE 


By:  R.  A.  Steinberg 
J.  J.  Rivera 


FLEET  SYSTEMS  DEPARTMENT 

THE  JOHNS  HOPKINS  UNIVERSITY  ■  APPLIED  PHYSICS  LABORATORY 

John*  Hopkkta  Road,  Laural,  Maryland  20707 
Opantlng  Under  Contract  N0Q03M7-C-S301  with  tha  Dapartmant  of  the  Navy 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


l 

'  \  t' 


''  {(" 


ABSTRACT 


A  mult/^resolution  signal  processing  algorithm  (*MRSI*)  was  devised  to  max¬ 
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Probabilities  of  detection,  determined  via  2500  image-based  Monte  Carlo 
simulations,  agree  with  well-known  analogous  results  of  psychophysical  ex¬ 
periments  performed  with  human  subjects.  A  new  method  is  presented  for 
supres/ing  false  alarms  caused  by  cloud  reflections  from  the  sea  surface. 
Regardless  of  whether  final  target  classification  is  performed  on-board  the 
missijfe  by  a  computer  algorithm  or  off-board  by  a  person  inspecting  imagery 
telemetered  from  the  missile,  MRSI  substantially  improves  predicted  system 
performance. 
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1.0  INTRODUCTION  AND  SUMMARY 


Self-guided  missiles  have  on  •  number  of  recent  oc¬ 
casions  demonstrated  high  effectiveness  in  operations 
against  surface  ships  [1].  Both  infrared  (IR)  and  radio 
frequency  (RF)  means  have  historically  each  found  ap¬ 
plication  as  the  sensing  element  in  antiship  missile 
(ASM)  guidance  units.  Some  well-known  advantages 
of  IR  technology  in  this  application  include  (2] 

•  passive  operation. 

•  good  resistance  to  jamming,  and 

•  high  spatial  resolution. 

The  impetus  for  high  spatial  resolution  follows  from 
potential  operational  needs  for  target  classification: 
assuming  adequate  signal  to-noise  ratio  (SNR),  high 
classification  accuracy  requires  high  spatial  resolution, 
regardless  of  whether  the  imagery  is  interpreted  by  a 
man [3 ,4]  or  processed  by  a  computer  [5].  Perhaps  the 
key  limitation  of  IR  ASM  seekers  is  limited  range  per¬ 
formance  under  conditions  of  poor  atmospheric 
visibility.  Thus,  the  purpose  of  the  signal  processor 
described  in  this  report  (referred  to  as  “MRSI”  for 
Multi-Resolution  Spatial  Integrator)  is  to  perform 
detection  of  targets  having  minimum  contrast  relative 
to  their  background,  in  a  digitized  two-dimensional 
image  The  MRSI  processing  approach  could  be  used 
to  optimize  the  detection  range  of  infrared  sensors 
against  lsuge  area  targets  such  as  ships. 

Means  previously  investigated  for  maximizing  IR 
sensor  SNR,  relevant  to  the  ASM  seeker  application, 
include  waveband  optimization  [6,7],  advanced  IR 
detector  developments  [8,9],  and  multi-frame  image 
processing  [10].  Nonetheless,  acquisition  range  remains 
an  important  performance  index,  and  improvements 
in  acquisition  range  an  important  objective,  in  the 
design  of  next  generation  antiship  seekers.  The  single- 
frame  signal  processing  approach  discussed  in  this 
report  is  complementary  to  earlier  approaches  for  SNR 
optimization. 

The  idea  for  our  new  signal  processing  concept  was 
anticipated  by  a  brief  study  indicating  that  detection 
ranges  obtained  by  human  observers  of  visual  displays 
could  under  a  broad  variety  of  conditions  greatly  ex¬ 
ceed  ranges  obtained  by  a  hot-spot  detection  algorithm 
[11].  (This  analysis  was  based  on  IR  ship  signatures 
cataloged  in  [12]  and  the  human  observer/display 
model  provided  in  [3,4].)  Thus,  the  predicted  perfor¬ 
mance  of  an  archetypal  “human  observer”  became  for 
us  the  standard  against  which  to  guage  the  perfor¬ 
mance  of  proposed  ship  detection  algorithms.  A 


similar  outlook  was  adopted  in  [13]  for  assessing  the 
performance  of  a  digitally-implemented  autotracker 
algorithm. 

Our  proposed  solution  to  the  ship  detection  problem 
is  comprised  of  a  separable  multi-resolution  niter, 
which  we  refer  to  as  a  Multi-Resolution  Spatial  In¬ 
tegrator  (MRSI),  first  described  in  [14]. 

An  early  multi-resolution  processing  application  in¬ 
volved  the  use  of  quad  trees  for  encoding  1-bit  images 
[IS].  A  variety  of  additional  multi-resolution  image 
processing  applications  and  computational  methods 
and  costs  are  discussed  in  [16,17].  The  multi-resolution 
processor  discussed  in  this  report  is  highly  specializ¬ 
ed  to  maximize  SNR  against  rectangular  and  nearly- 
rectangular  targets  of  known  orientation,  but  whose 
size,  aspect  ratio,  and  position  in  the  image  are 
unknown.  The  MRSI  approach  is  by  design  ideally 
suited  for  processing  data  from  a  linear  array  scann¬ 
ing  sensor,  i.e.,  image  data  acquired  sequentially  by 
column. 

The  time  domain  responses  and  SNR  enhancement 
provided  by  MRSI  processing  are  amenable  to  exact 
analytical  computation  for  only  a  few  geometrically 
simple  target  shapes.  A  number  of  such  analytically 
predictable  processor  responses  was  used  to  validate 
a  general  image-based  computer  simulation.  Tbst 
targets  used  for  this  purpose  included  point  targets, 
bar  targets,  white  gaussian  noise,  and  bar  targets  in 
noise.  The  computer  simulation  was  then  exercised 
against  a  variety  of  ship  images  representative  prin¬ 
cipally  of  frigates  and  aircraft  carriers,  seen  from 
ranges  of  from  S  to  30  nautical  miles  (nmi)  at  viewing 
aspects  from  90°  (beam)  to  10 Threshold  SNRs  for 
detection  were  established  in  simulation  by  adding  to 
each  test  image  progressively  increasing  levels  of  noise. 

It  was  found  that  MRSI  processing,  under  a  broa  ' 
variety  of  conditions,  provides  about  ten-fold  SNk 
enhancement  (“processing  gain”)  relative  to  hot-spot 
detection.  Moreover,  the  processing  gain  turns  out  to 
have  been  computable,  with  no  more  than  about  104k 
maximum  error  in  every  c'se,  as 

G  =  riN?lx  (1) 

where  NPIX  is  the  number  of  pixels  occupied  by  the 
target  in  the  image,  and  »j  is  given  by 

n  =  0.81.  (2) 

The  processing  constant,  t>,  is  remarkably  invariant 
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with  respect  to  sensor/ship  range  (i.e.,  scale),  viewing 
aspect,  and  ship  type. 

The  existance  of  a  simply-computed  analytical  ex¬ 
pression  for  processing  gain,  Eq.  1,  is  in  itself  a  signifi¬ 
cant  result,  as  it  eliminates  the  need  to  perform  image- 
based  simulations  in  connection  with  cost/perfor¬ 
mance  trade  studies  involving  such  additional  aspects 
of  design  as  focal  plane  sensitivity  requirements.  Also, 
an  expression  for  ij  has  been  derived  that  can  be 
evaluated  simply  from  ship  profiles,  such  as  those  pro¬ 
vided  in  Janets  Fighting  Ships  [18].  Consequently,  our 
results  are  readily  generalized  to  ship  types  other  than 
those  for  which  detailed  simulations  have  been  con¬ 
ducted.  At  least  as  regards  first-order  system  perfor¬ 
mance  assessments,  Eq.  1  may  be  used  as  a  vastly 
simpler  substitute  for  detailed  image-based 
simulations. 

The  principal  simulation  results  are  plots  of  detec¬ 
tion  probability  (P„)  vs.  SNR  (cf.  Figs.  42  and  43), 
based  on  2500  image-based  Monte  Carlo  trials.  A 
curve  fit  to  the  numerically-developed  probabilities  of 
detection  appears  to  provide  a  universal  curve  that  can 
be  used  predictively  for  ships  of  differing  class,  observ¬ 
ed  at  a  variety  of  aspects  and  distances.  Interestingly, 
the  curve  fit  to  MRSI/Monte  Carlo  PD  results  coin¬ 


cides  with  a  previously  published  curve  fit  to  data  ob¬ 
tained  from  visual  detection  experiments  performed 
with  human  observers  of  TV  displays  (3). 

TWo  ASM  system  concepts  are  discussed  in  Section 
15.0,  with  each  concept  employing  a  different  means 
of  target  classification,  viz.,  classification  performed 
on-board  the  missile  by  an  autonomous  target  recogni¬ 
tion  (ATR)  algorithm  or  off-board  by  a  person  inspec¬ 
ting  imagery  telemetered  from  the  missile.  Using  MRSI 
to  cue  a  slow  scan/high  sensitivity  "classification 
mode",  as  proposed  by  W.  J.  TYopf  (APL),  results  in 
substantial  improvements  in  predicted  performance  for 
both  system  concepts. 

In  Section  16.0  we  present  a  new  method  for  sup¬ 
pressing  false  alarms  caused  by  cloud  reflections  from 
the  sea  surface. 

The  image-based  simulation  itself  is  coded  in  FOR¬ 
TRAN  and  operates  on  ship  images  of  64  x  256  pix¬ 
els,  at  35  spatial  resolutions  (App.  A).  The  program 
runs  both  on  the  APL  National  Advanced  Systems 
(NAS)  mainframe  and  on  an  IBM  PC  AT  computer, 
and  requires  about  450K  bytes  of  memory.  All  simula¬ 
tion  results  provided  in  this  report  (more  than  2800 
processed  images)  were  developed  on  an  IBM  PC  AT. 
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2.0  ASSUMPTIONS  AND  DEFINITIONS 


Scan 

Olfaction 


Figure  1.  Scanning  IR  Makar  concapt.  Vortical  linear 
array  of  I  infrared  detectors  la  scanned  In  azimuth, 
generating  I  digital  waveforms  which  are  then  Input  to 
the  MRSI  processor.  M RSI's  output  is  a  list  of  objects, 
each  characterized  in  terms  of  Its  position,  dimensions 
(az  x  el)  and  brightness. 


it  is  assumed  that  the  MRSI  signal  processor  is  at¬ 
tached  to  an  imaging  sensor  that  performs  surveillance 
of  a  scene  once  very  TF  seconds,  where  the  interval 
Tr  is  called  the  sensor's  "frame  time?*  (Fig.  1). 

The  angular  dimensions  of  a  single  detector's  field- 
of-view,  denoted  as  a  (in-scan)  ano  0  (cross-scan)  are 
collectively  referred  to  as  the  seeker’s  “spatial  resolu¬ 
tion”  or  “instantaneous  f»eld-of-view"  (IFOV)  and  are 
specified  in  milliradian  (mr)  units.  The  voltage,  V,(j)> 
read  out  from  detector  i  at  sampling  time  j  is 
presumably  proportional  to  the  integral  of  the  scene 
radiance  over  a  rectangular  region  of  angular  dimen¬ 
sions  a  X  0. 

The  viewed  scene  is  assumed  to  be  a  rectangle  of 
angular  dimensions  A  x  B,  mr.  The  sensor  IFOV  sub¬ 
divides  the  total  scene  into  a  rectangular  grid  of 
samples,  each  of  angular  dimensions  a  x  &  mr.  The 
total  scene  (also  called  a  “frame”)  and  the  sampling 
geometry  are  depicted  in  Fig.  2.  Also  shown  in  Fig. 
2  is  a  rectangular  coordinate  frame  (s,x).  Samnles  are 
obtained  on  a  uniform  rectangular  grid,  with  adjacent 
samples  separated  by  angle  a/u  in  the  s-direction  and 
by  angle  0/v  in  the  x-direction.  (0/v  is,  of  course,  the 
angular  separation  of  adjacent  detector  elements  in 
the  detector  array.)  As  mentioned  ibove,  each  sample 


is  considered  to  be  an  integral  of  the  scene  brightness 
over  a  small  rectangular  region  of  angular  dimensions 
a  x  A  As  shown  in  Fig.  2,  dimensions  (A,  a,  a/u) 
are  parallel  to  the  s-axis,  while  dimensions  (A,  d,  0/v) 
are  parallel  to  the  x-axis. 

Depending  on  how  the  sensor  is  implemented,  the 
spatial  samples  may  be  gapped,  overlapped,  or  precise¬ 
ly  adjacent.  Figure  2  depicts  an  example  for  which  the 
samples  are  gapped  in  both  the  a-  and  x-directions, 
ije.,  for  which 


u,  u  <  1.  (3) 

Samples  are  precisely  adjacent  in  the  s-direction  when 

u  *  1.  (4) 

Samples  are  overlapped  in  the  s-direction  when 

u  >  l.  (5) 

Analogous  to  Eqs.  4  and  S,  samples  are  precisely  ad¬ 
jacent  and  overlapped  in  the  x-direction  when  v=! 
and  u>l,  respectively. 

In  general,  the  ratios  of  sample  size  to  sample  spac¬ 
ing,  u  and  v,  are  different  in  value.  For  example;  depen- 
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Figure  2.  Total  scene  la  s  rectangle  of  angular  dimen¬ 
sions,  AxB.  Spatial  samples  are  of  angular  dimensions, 
a  x  0.  Samples  are  obtained  on  centers  displaced  by 
angles  a/u  and  0/v.  All  angles  are  measured  In  mllllra- 
dlana  (mr.) 
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ding  on  the  sensor’s  implementation,  samples  may  be 
adjacent  in  s  and  gapped  in  x(u  - 1 ,  v  <  1 )  or  overlap¬ 
ped  in  s  and  adjacent  in  x(u>  1,  u=l);  etc. 

In  the  following  discussion,  a  scanning  sensor  im¬ 
plementation  is  assumed,  in  winch  the  s-coordinate 
(Fig.  2)  coincides  with  the  direction  of  scan,  and  the 
x-coordinate  is  variously  referred  to  as  the  “cross-scan” 
or  “cross-channel”  coordinate.  Quantity  u  is  conven¬ 
tionally  referred  to  as  the  in-scan  sample  factor,  hav¬ 
ing  units  of  samples  per  dwell. 

Although  the  following  discussion  adopts  the 
parlance  of  scanning  sensors,  the  MRSI  processor  is 
readily  adapted  to  staring  mosaic  sr  .isors  as  well. 

Regardless  of  the  numerical  values  of  the  parameters 
defined  in  Fig.  2,  the  sensor  develops  every  TF 
seconds  a  matrix  of  numbers,  (  Vt(J)  I,  in  ’vhich  each 
matrix  element  is  proportional  to  the  brightness  of  the 
scene,  averaged  over  a  neighborhood  of  dimensions 
a  x  0,  centered  on  the  angular  coordinates  (Ja/u , 
i0/v). 

The  detector  index,  i,  is  also  the  cross-channel  im¬ 
age  coordinate,  taking  on  the  range  of  values 

1  sis/  =  [Bv/0],  (6) 

where  square  bracktts  (•]  denote  the  integer  part  of 
their  argument.  Similarly,  the  discrete  time  variable, 
j,  takes  on  the  range  of  values 


Figure  3.  Part  a):  The  field-of-view  of  an  IR  detector 
is  scanned  across  an  ideal  bar  target.  T  denotes 
temperature.  Detector  and  target  dimensions,  a  and  y 
respectively,  are  in  angle  units  (mrad).  Angular  rate  of 
scan  is  a  constant,  0,  mrad/sec. 

Part  b):  Appearance  of  the  voltage  waveform  obtained 
at  the  detector's  output  for  the  situation  depicted  In 
part  a). 


1  sj*J=  [Au/aj.  (7) 

During  an>  given  sampling  interval  (i.e.,  for  a  given 
value  of  J),  the  MRSI  processor’s  input  is  a  vector  Vi 
of  dimension  /.  Equivalently,  the  processor  operates 
upon  an  input  consisting  of  /  digital  waveforms. 

If  the  IR  seeker  were  positioned  in  the  far  field  of 
a  bar  target,  the  voltage  output  of  each  detector  chan¬ 
nel  would  appear  approximately  as  shown  in  Fig.  3. 

For  the  purposes  of  this  discussion,  an  object  of 
angular  dimensions: 

ys,  in  the  direction  of  scan  (“in-scan”) 

yx,  cross  scan, 

is  said  to  be  “resolved  in-scan”  when  ys>ot,  and 
“resolved  cross-scan”  when  yx>0.  An  objec'.  is  said 
to  be  “fully  resolved”  when  it  is  resolved  both  in-scan 
and  cross-scan. 

As  discussed  in  [11],  ship  targets  arc  fully  resolved 
at  initial  detection,  for  all  cases  of  practical  interest. 

The  “channel  signal-to-noise  ratio”,  SNRC,  is 
calculated  as  (cf.  Fig.  3b) 

SNRC  =  ( Vh  -  Vc)/Vrm  (8) 


i.e.,  as  the  ratio  of  peak  signal  to  root-mean-square 
(RMS)  noise,  at  the  output  of  a  single  detector  channel. 

A  basic  measure  of  seeker  sensitivity  is  its  “noise 
equivalent  temperature”,  NEAT,  calculated  as  (cf. 
Fig.  3) 

NEAT  =  (Th  -  Tc)/SNRC  (9) 

with  SNRC  given  by  Eq.  8. 

The  sensitivity  parameter  NEAT  may  be  calculated 
using  measurements  made  against  calibrated 
laboratory  targets,  using  Eqs.  8  and  9.  Quantity  NEAT 
may  also  be  calculated  (or  predicted  prior  to  sensor 
fabrication)  as  a  function  of  such  key  seeker  optical 
characteristics  as  aperture,  $,  a,  detector  detectivity 
(D*),  etc,  [19]. 

The  number  of  pixels  (spatial  resolution  elements) 
occupied  by  a  rectangular  target  of  size  (in-scan) 
x  7v  (cross-scan)}  is  given  by 

npix  =  0°) 

Equation  9  is  presented  above  from  the  perspective 
of  using  calibrated  bar  target  data  to  obtain  a  measured 
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va'ue  of  NEAT,  as  a  means  of  characterizing  seeker 
sensitivity.  Once  NEAT  has  been  obtained,  a  variant 
of  Eq.  9  may  then  be  used  to  predict  the  SNR  expected 
in  a  tactical  scenario,  viz., 

SNRC  =  AT/NEAT  01) 

where  the  quantity  A  T  is  now  the  “apparent”  contrast 


temperature  of  the  tactical  target.  Apparent  ship  con¬ 
trast  temperature,  A  T,  is  sensitive  to  a  wide  variety  of 
parameters  characteristic  of  the  target,  the  background, 
and  the  atmosphere  (20]. 

Quantities  AT  and  SNRC  are  most  generally  deter¬ 
mined  as  spatial  distributions  since  they  may  take  on 
different  values  for  each  pixel  in  the  frame. 
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3.0  HUMAN  VISION  AS  A  MODEL  OF  PERFECT  PROCESSING 


At  least  as  regards  spatial  resolution  (a,  /S),  re¬ 
quirements  for  target  recognition  are  at  odds  with  re¬ 
quirements  for  long  range  inquisition.  Assuming  ade¬ 
quate  SNR,  high  accuracy  in  target  recognition  requires 
high  spatial  resolution,  regardless  of  whether  the  im¬ 
agery  is  interpreted  by  a  man  [3,4]  or  processed  by  a 
computer  [5].  Thus,  the  objectives  of  recognition  are 
furthered  by  making  the  number  of  pixels  on  target 
(Nplx)  very  large.  On  the  other  hand,  SNR  is  max¬ 
imized  when  the  resolution  element  is  matched  in  size 
to  the  target,  i.e.,  acquisition  range  is  maximized  when 
N„x  =  1.  . 

The  apparent  discrepancy  in  resolution  requirements 
for  IR  acquisition  and  recognition  can  be  resolved  via 
detection  processing  modeled  after  our  understanding 
of  eye/brain  processing  in  human  vision  [3]. 

In  performing  detection  of  resolved  targets  in  elec¬ 
tronic  images,  humans  are  apparently  capable  of  per¬ 
forming  adaptive  integration  over  the  target’s  area  - 
an  astounding  feat  in  light  of  the  fact  that  ‘he  boun¬ 
daries  for  area  integration  must  be  established  prior 
to  conscious  detection. 

The  remarkable  human  capacity  for  subliminal  in¬ 
tegration  over  arbitrarily  sized  and  located  target  areas 
is  referred  to  by  Rosell  as  “the  perfect  synchronous 
integrator  model”  of  human  vision,  and  is  attributed 
by  Rosell  to  original  discovery  by  Otto  Schade,  Sr.  [4]. 

The  effect  of  perfect  synchronous  integration  is  to 
enhance  by  N”IX  the  single  channel  SNR,  SNR,-, 
given  by  Eq.  11.  A  further  characteristic  of  human  vi¬ 
sion  is  that  the  process  of  detection  cannot  be 
separated  from  the  process  of  estimating  target  angular 
dimensions;  the  two  processes  are  performed  in  a 
coupled  and  purely  parallel  fashion. 

Figure  4  depicts  a  conceptual  model  of  the  human 
vision  system  (II VS)  detection  process  based  on  a 
parallel  array  of  spatial  matched  filters,  in  which  each 
filter  in  the  array  corresponds  to  a  possible  target 
shape.  Assuming  that  the  HVS  is  indeed  a  perfect  syn¬ 
chronous  integrator  (PSI),  every  possible  target  shape 
is  represented  in  the  parallel  filter  bank,  as  well  as  all 


Figure  4.  Conceptual  model  of  the  human  vision 
system  (HVS)  detection  process,  based  on  an  infinite¬ 
dimensional  array  of  spatial  matched  filters. 

variants  of  each  shape  obtainable  by  the  processes  of 
translation,  rotation,  and  scaling. 

The  detection  model  diagrammed  in  Fig.  4  is  not 
directly  amenable  to  simulation  or  digital  realization, 
since  there  are  an  infinite  number  of  filters  in  the  ar¬ 
ray.  We  nonetheless  fearlessly  assume  that  it  is  possi¬ 
ble  to  implement  with  digital  hardware  an  approxima¬ 
tion  to  synchronous  area  integration,  which  leads  to 
the  definition  of  a  detection  SNR, 

SNR,',  =  H  •  SNRC  •  N»x  .  (12) 

Equivalently,  the  effect  of  synchronous  area  integra¬ 
tion  can  be  expressed  in  terms  of  a  “processing  gain” 

G  =  (SNRdel/SNRc)  =  t,  Nplx  ,  (13) 

presented  previously  as  Eq.  1.  Quantity  t/  :n  Eqs.  12 
and  13  is  a  “processor  efficiency”  that  is  included  to 
allow  for  the  possibility  of  suboptimal  processor  per¬ 
formance.  In  the  usage  of  this  report,  a  PSI  is  any  pro¬ 
cessor,  human  or  machine,  characterized  by  the  pro¬ 
perty  rj  =  1.  By  extension,  an  imperfect  synchronous 
integrator  (ISI)  is  a  processor  for  which  r?  <  1. 
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4.0  ACQUISITION  RANGE  CALCULATIONS 


It  follows  from  Eqs.  11  and  12  that 

AT  =  r'  ■  SNRdtl  ■  NEAT  /  Nplx  .  (14) 

Assuming  somewhat  arbitrarily  a  threshold  value  for 
SNRdet 


where 

AShiP  =  ship’s  projected  area,  mJ 

a,  P  =  sensor  pixel  dimensions  (cf.  Fig.  2),  mr 

R  -  range  from  seeker  to  ship,  km. 


SNRTdtl  =  5  ,  (15) 

the  value  of  AT  required  for  detection,  RAT 

RAT  =  5  •  NEAT / i)Np,x  (16) 


is  obtained. 

Quantity  RAF  is  similar  to  the  Minimum  Resolvable 
Temperature  (MRT)  parameter  conventionally  used  to 
describe  the  performance  of  human  observers  of  ther¬ 
mal  image  displays.  In  fact,  RAT  can  be  made  equal 
to  MRT  if 


•  SNRTdt ,  is  set  equal  to  2.8  (rather  than  5),  cor¬ 
responding  to  the  experimentally  determined 
value  required  by  human  observers  to  achieve 
50%  probability  of  detection, 

•  r)  is  set  equal  to  unity,  and 

•  an  additional  factor  (TF/0.2)''  is  included, 
where  TF  is  frame  time  (typically,  1/30  se¬ 
cond),  to  account  for  the  smoothing  of  rapid 
temporal  variations  by  the  eye’s  limited  response 
time  (~0.2  s). 


The  application  of  MRT-like  descriptors  to 
automatic  target  screeners  is  also  proposed  in  [21]. 

In  the  absence  of  spatial  integration  processing  (i.e., 
performing  threshold  detection  directly  on  the  in¬ 
dividual  channel  outputs),  Eq.  16  becomes 

RAT  -  5  •  NEAT  .  (17) 


The  ratio  of  Eqs.  16  and  17  is,  of  course,  Eq.  13. 

For  fixed  sensor  IFOV  (a  and  P  in  Fig.  2),  NPIX 
decreases  with  increasing  range.  Hence,  from  Eq.  16, 
greater  apparent  contrast  is  required  to  detect  distant 
targets  than  is  required  to  detect  closer  targets.  Ap¬ 
parent  contrast  required  for  simple  hot  spot  detection, 
given  by  Eq.  17,  is  independent  of  NP1X  and  thus  is 
also  independent  of  target  range. 

Tb  further  evaluate  Eq.  16,  the  number  of  pixels  on 
target,  NP1X,  is  expressed  as 


^p/x  Ash ip/ct0R2 


From  Eqs.  13  and  18 

G  =  r,A«ip/R{a0 )*  .  (19) 

The  apparent  contrast  required  for  detection  may  be 
written  in  terms  of  Eq.  19  as 

RAT  =  5  •  NEAT/G  .  (20) 


Equations  17,  19,  and  20  can  be  used  for  calculating 
seeker  acquisition  range  against  ship  targets,  as  in¬ 
dicated  graphically  in  Fig.  5,  for  the  strawman  seeker 
and  target  characteristics  provided  in  Thble  1. 

Curve  A  in  Fig.  5  is  RAT  for  a  simple  hot  spot  detec¬ 
tion  algorithm,  from  Eq.  17.  Curves  B  and  C  are  RAT 
for  synchronous  integration  processing,  from  Eq.  20, 
for  the  front-  and  beam-aspect  targets,  respectively. 

Also  shown  in  Fig.  5  (curve  D)  is  a  representative 
ship  “signature”  (i.e.,  ship  apparent  contrast 


Range  to  ship 


Figure  5,  Acquisition  range  calculations  for 
strawman  seeker  and  target  characteristics  provided 
In  Table  1. 


(18) 
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TabU  1 

Strawman  saeker  and  target  parameters 


Seeker 

Parameter 

Meaning 

Assumed 

Value 

a 

In-scan  resolution 

0.2  mr 

p 

Cross-scan  resolution 

0  2  mr 

V 

Processor  efficiency 

1 

NEAT 

Temperature  sensitivity 

0.1  K 

Thrget 

Parameter 

Meaning 

Assumed 

Value 

Ashlp  (90°)  Ship  beam-aspect  area  4750  m: 

Aihlp  (0°)  Ship  bow-aspect  area  1100  mJ 


temperature  vs.  range)  for  a  particular  locale,  season, 
time  of  day,  probability  of  weather,  sensor  height- 
above-water,  and  IR  waveband  [12]. 

Intersection  points  /,  and  /2  in  Fig.  5  indicate  that 
synchronous  integration  processing  achieves  acquisi¬ 
tion  ranges  of  7.5  and  10.3  kilometers  (km)  against 
the  front-  and  beam-aspect  targets,  respectively.  Lack 
of  intersection  between  curves  D  and  A  indicates  that 
hot  spot  detection  fails  to  achieve  target  acquisition. 

At  this  juncture,  it  should  be  obvious  that  syn¬ 
chronous  area  integration  offers  considerable  perfor¬ 
mance  advantage  relative  to  hot  spot  detection.  Less 
obvious  is  how  these  benefits  are  to  be  obtained,  since 
the  only  model  presented  thus  far  for  a  synchronous 
area  integrator  (viz.,  Fig.  4)  is  unrealizable. 


* 


* 


1 

f 
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5.0  AN  IMPERFECT  SYNCHRONOUS  INTEGRATOR 


In  this  section  we  develop  the  structure  of  a 
realizable  signal  processor  (an  “imperfect  synchronous 
integrator’  or  ISI)  having  performance  characteristics 
almost  as  good  as  the  unrealizable  PSI  depicted  in  Fig. 
4. 

We  start  the  ISI  design  with  the  structure  of  Fig. 
4,  i.e.,  a  processor  configured  as  a  parallel  bank  of 
spatial  filters.  The  PSI  is  able  to  accommodate  target 
profile  uncertainties  in  shape,  scale,  position,  and 
orientation.  The  first  step  in  reducing  processor  com¬ 
plexity  is  to  note  that  the  rotational  degree  of  freedom 
is  unneeded  in  the  present  application:  to  sea-skimming 
missiles  flying  wings-level,  the  major  axis  of  ship 
targets  will  always  appear  parallel  to  the  horizon. 

The  next  step  in  reducing  processor  complexity  is 
to  constrain  the  spatial  filter  responses  to  be  rectangles. 
Thus,  when  the  original  image  consists  of  a  single 
bright  pixel  against  an  otherwise  uniform  background, 
the  output  image  of  each  filter  will  be  a  rectangle  of 
dimensions  k  pixels  (cross-scan)  by  n  pixels  (in-scan). 
Consequently,  each  filter  in  the  filter  bank  can  be  label¬ 
ed  with  its  “resolution  indices”,  k  and  rt,  and  general¬ 
ly  described  in  terms  of  “the  (k,n)  filter!’  The  resulting 
suboptimal  approximation  of  Fig.  4  is  depicted  in 
Fig.  6. 

Finally,  the  required  number  of  filters  in  the  filter 
bank  is  further  limited  by  constraining  the  resolution 
indices  to  be  powers  of  2,  between  specified  minimum 
and  maximum  values  based  on  the  likely  target  angular 
dimensions  at  initial  detection.  For  example,  one  may 
choose 


corresponding  to  a  filter  bank  comprised  of  15  spatial 
filters.  The  number  of  spatial  filters  may  be  reduced 
still  further  by  culling  out  filters  that  correspond  to 
unlikely  target  aspect  ratios.  As  indicated  in  Fig.  7  for 
example,  there  may  be  very  little  penalty  in  removing 
filters  corresponding  to  (k,n)  =  (8,4)  and  (k,n)  = 
(2,64);  targets  having  such  aspect  ratios  can  still  be 
detected,  although  with  reduced  SNR  gain. 

The  parallel  bank  of  13  filters  corresponding  to  Fig. 
7  may  be  conceived  as  having  a  single  image  as  input 
and  generating  13  blurred  /reduced  resolution  images 
as  its  output.  (In  the  image-based  simulations  discuss¬ 
ed  in  Section  14,  a  processor  having  35  spatial  filters 
was  routinely  simulated  on  an  IBM  PC  AT.)  It  is  for¬ 
tunate  from  the  standpoint  of  data  throughput  re¬ 
quirements  that  the  sampling  rate  for  each  of  the  13 
reduced  resolution  images  is  substantially  lower  than 
for  the  original  full  resolution  image. 

Considering  Fig.  6,  n  —  4  is  the  lowest  resolution 
at  which  the  scene  need  be  observed  from  the  stand¬ 
point  of  matched  filtering  for  detection.  This  has  sug¬ 
gested  the  structure  of  a  detection  preprocessor 
depicted  in  Fig.  8.  Each  of  the  /  detector  channels  is 
processed  by  an  identical  such  preprocessor.  The  nota¬ 
tion  for  image  data  generally  is 

K"U)  (22) 

where  we  recall  that  indices  (ij)  refer  to  position  in 
the  image  (cf.  Fig.  1)  and  indices  (k,n)  denote  the 
spatial  resolution  of  the  image.  Thus 


v';'U) 

in  Fig.  8  is  the  image  data  as  originally  output  from 
the  detector  array,  at  the  original  (k,n)  =  (1,1) 


Figure  8.  Imperfect  synchronous  Integrator  (ISI)  ob-  Figure  7.  Matrix  of  spatial  resolutions  for  a  spatial 
talned  from  Figure  4  by  selecting  spatial  filter  filter  array, 
responses  as  unrotated  rectangles. 
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Figure  8.  Detection  preprocessor  block  diagram. 

resolution.  As  shown  in  Fig.  8,  the  voltage  time 
waveform  from  detector  i,  denoted 

v‘JU) 

(where  j  is  time)  is  smoothed  by  a  boxcar  integrator 
and  decimated  four-fold  in  time;  i.e.,  only  one  of  each 
four  samples  output  from  the  boxcar  integrator  is  us¬ 
ed  in  subsequent  signal  processing  stages.  The  effect 
of  the  detection  preprocessing  step  is  to  provide  an  im¬ 
mediate  two-fold  improvement  in  SNR  and  an  im¬ 
mediate  four-fold  reduction  in  subsequent  data 
throughput  requirements. 

Of  course,  t’.ie  two- fold  SNR  improvement  is  only 
experienced  for  targets  at  least  four  samples  in  in-scan 
extent;  smaller  targets  will  have  less  SNR  improvement, 
with  a  single-sample  target  actually  experiencing  a  two¬ 
fold  degradation  in  SNR.  This  drawback  is  more  ap¬ 
parent  than  real  for  the  following  two  reasons. 

•  Subsequent  to  detection,  it  is  necessary  to 
transfer  processing  to  a  classification  step. 
Targets  occupying  s4  samples  in-scan  are  too 
small  to  permit  image-based  classification, 
e  Size  is  a  key  indicator  of  ship  target  value.  Small 
ships  are  generally  of  less  interest  as  targets  than 
are  large  ships. 

The  input /output  relationship  for  the  Fig.  8  pre¬ 
processor  may  be  expressed  as 

V)AU)  =  '/<  tQ  V]'(4j-k)  ,  (23) 

j— 1,2,...  [Au/4ot]  . 

The  maximum  value  of  j  in  Eq.  23 

j  max  =  lAu/4a]  (24) 

is  only  one-fourth  as  large  as  jmax  in  Eq.  7  due  to  the 
four-fold  decimation  step. 

Writing 

A  =  QTf  (25) 

where 


A  =  total  in-scan  field-of-view,  mr 
0  =  angular  rate-of-scan  mr/s 
7>  =  frame  time,  s 

and  recognizing  that  the  sample  rate  per  channel,  fs 
(Hz),  is  given  by 

=  JmJTF  (26) 

from  Eqs.  (24)-(26),  it  can  be  found  that 

fs  =  GuMa  (27) 

The  factor  4  in  the  denominator  of  Eq.  27  is  due  to 
the  decimation  step  in  Fig.  8. 

The  factor  (1/4)  preceding  the  sum  in  Eq.  23  is  a 
normalization  that  may  not  be  required  in  a  hardware 
realization  of  the  processor. 

The  general  structure  of  the  complete  ISI,  which  is 
called  a  Multi-Resolution  Spatial  Integrator  (MRSI), 
is  shown  in  Fig.  9.  The  elements  of  MRSI  are: 

•  the  detection  preprocessor  (Fig.  8), 

•  cross-channel  averaging  (XAV,  Section  6), 


Detector  Detection  Filter  band 


Detections 


Figure  9.  MRSI  signal  processor  block  diagram. 
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•  in-scan  averaging  (IAV,  Section  7), 

•  background  estimation  and  subtraction  (BES, 
Section  9),  and 

•  noise  estimation  and  threshold  /decision  (NE 
and  T/D,  Section  10). 

Each  capitalized  acronym  (XAV,  IAV,  etc.)  is  the  name 
of  the  subroutine  that  implements  the  corresponding 
processor  function  in  the  FORTRAN  simulation 
(App.  A). 

The  action  of  the  detection  preprocessor  has  already 
been  discussed  in  connection  with  Eq.  23  and  Fig.  8. 
The  other  processing  blocks  in  Fig.  9  will  be  discuss¬ 
ed  in  subsequent  sections  of  this  report. 

Although  the  detector  array  depicted  in  Fig.  9  is 
comprised  of  square  non-staggered  elements  with  a 
small  cross-channel  gap,  the  MRS1  processor  is  readily 
applied  to  arrays  having  all  possible  combinations  of 
the  following  attributes:  arbitrary  values  of  v;  stag¬ 
ger  offset  between  adjacent  detector  channels;  and 
time-delay-and-integration. 

The  means  for  implementing  detection  preprocess¬ 
ing  as  in  Fig.  8  is  an  array  of  1  identical  processing 
modules  operating  simultaneously  and  uncoupled;  i.e., 


in  a  purely  parallel  fashion.  A  means  for  implemen¬ 
ting  an  array  of  spatial  filters  at  varying  spatial  resolu¬ 
tions  is  described  in  Sections  6  through  8.  Like  the 
preprocessor,  Fig.  8,  the  multi-resolution  processor 
comprises  a  large  number  of  simple  and  identical  pro¬ 
cessing  elements  operating  in  parallel.  The  spatial  filter 
algorithm  is  separable  in  in-scan  and  cross-scan  opera¬ 
tions,  with  the  cross-scan  operations  being  perform¬ 
ed  first.  Thus,  we  operate  on  the  detection  preprocessor 
outputs,  viz., 

V)'*U)  . 

with  a  “cross-channel  averaging  processor”,  XAV,  to 
obtain 

K*U)  .  *  =  2,4,8 

which  are  then  provided  as  input  to  an  "in-scan  averag¬ 
ing  processor”,  IAV,  to  ob  ain  the  complete  set  of 
reduced  resolution  images, 

V*'"U)  ,  k  =  2,4,8 

n  =  4,8,16,32,64 
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6.0  CROSSCHANNEL  AVERAGING  PROCESSOR  (XAV) 


The  structure  of  the  cross-channel  averaging  pro¬ 
cessor,  XAV,  is  shown  in  Fig.  10.  Although  three  stages 
of  cross-channel  processing  are  required  to  obtain  the 
three  levels  of  cross-channel  resolution  indicated  in  Fig. 
7  ( k  —  2, 4,  and  8),  for  simplicity  of  illustration,  only 
two  of  the  stages  are  depicted.  Again  for  simplicity, 
in  Fig.  10,  the  four  quantities  A,  B,C,  and  D  are  defin¬ 
ed  as  follows: 

A  =  (28) 

B  =  V'SU) 

C  =  V\*(j) 

D  =  V\'*U) 

Through  the  first  stage  of  cross-channel  averaging,  we 
obtain 

V]A  U)  =  (A  +  B)/2  =  VilK',*4  U)  +  K*-4  C/)  1  (29) 

viA  u)  =  (B+o/2  =  'AtFi'4o')+ v;-4o)] 

VY  O')  =  (C+D)/ 2  =  Vi l Kj'4  (/)  +  »^i'4  U)  1 

Generalizing  the  relationship  we  see  developing  in  Eq. 
29, 

VYU)  =  >/i  V};4r.,0)  ,  (30) 

f  — 1,2,...  (/-I) 


descriptive  of  the  action  of  the  first  XAV  stage. 

Following  the  XAV  second  stage  logic  in  Fig.  10  we 
obtain, 

4 

V*AU)  =  (A  +  B+C+D)/ 4  *  '/4  E  VYU)  . 
(=1,2,...  (7-3) 

which  generalizes  to 

K'U)  =  i  V'Yr-lU)  ,  (31) 

i=l,2,...  (7-3) 

The  output  of  the  XAV  third  stage  (not  shown  in 
Fig.  10)  is  found  to  be 

V*Y(j)  =  1/8  E  V)Ar_ ,  (!)  ,  lJ2) 
i— 1,2,.. '  \l  ' )  . 

Inspection  of  F;qs  ?G  through  32  indicates  the 
general  input/output  description  of  XAV  can  be  writ¬ 
ten  as 

►TO)  =  1  /kEx  F;;4r.,0)  .  (33a) 

A:  =  2,4,8 

(  =  1,2,...  (7 — A -4- 1) 


Figure  10.  Cross-channel  averaging  processor,  XAV.  Inputs  to  XAV,  denoted  Vi1  >*()),  are  provided  by  detection 
preprocessor  (Figure  8).  Outputs  of  XAV,  denoted  by  unterminated  arrows  pointing  downward,  are  provided  to 
IAV  for  further  processing. 
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Cross-channel  averaging  successively  reduces  the 
number  of  elevation  channels  (and  hence  the  max¬ 
imum  value  of  0  at  progressively  higher  resolutions 

/(*)  m  max(i }  =  /-*  +  l  ,  (34) 

where  the  number  of  elevation  channels  at  the  original 
resolution,  /,  is  given  by  Eq.  6. 

It  follows  from  Eq.  33  that  cross-channel  process¬ 
ing  is  instantaneous,  i.e.,  determination  of  V*  A  at 
time  j  depends  only  on  values  of  K^4 ,  also  at  sampl¬ 
ing  time  j. 


The  position  of  XAV  within  the  overall  processing 
structure  is  indicated  in  Fig.  9. 

Equivalent  to  Eq.  33a,  we  can  write 

V?AU)  =  'AlF*-40‘)  +  KUU))  ,  (33b) 
*  =  1,2,4 

i=l,2,...  (/— 2*+  1)  . 

The  FORTRAN  routine  XAV  (App.  A)  implements 
Eq.  33b. 
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7.0  IN-SCAN  AVERAGING  PROCESSOR  (IAV) 


The  in-scan  averaging  processor,  (IAV),  is  comprised 
of  a  large  number  of  identical  processing  “sections", 
each  structured  as  shown  in  Fig.  11.  Since  a  separate 
IAV  section  is  required  for  each  distinct  pair  of  values 
(i,k),  the  total  number  of  required  sections  (each  iden¬ 
tical  to  Fig.  11)  may  be  computed  as 

n,av  =  E  Hk) 

With  I{k)  given  by  Eq.  34, 

Niav  =  E  (I—k+ 1 )  =  3/- 1 1  .  (35) 

**2,4,8 

For  example,  if  the  scanning  seeker  has 

/  =  128  (36) 

detector  channels,  from  Eq.  35,  it  is  found  that 

Niav  =  373  (37) 

IAV  sections,  each  identical  to  Fig.  11,  will  be  required. 

More  generally,  the  number  of  required  IAV  sections 
may  be  simply  estimated  as  /  times  the  number  of  re¬ 
quired  cross-scan  averaging  stages. 

The  action  of  IAV  is  developed  with  reference  to  Fig. 
11  as 


V,M(|)  -  from  XAV 


Figure  11.  Transversal  filter  realization  of  one  section 
of  the  in-scan  averaging  processor,  IAV.  Inputs  to  IAV 
are  provided  by  the  cross-scan  averager  (XAV).  A  large 
number  of  such  sections,  operating  simultaneously 
and  uncoupled,  comprise  the  complete  IAV.  Output  of 
IAV  is  an  array  of  reduced  resolution  images.  Process¬ 
ing  elements  denoted  as  z-1  are  unit  delays.  Final  8 
delay  elements  and  final  summing  node  are  not  shown. 


2 

V*'\j)  =  Zi  L  V*4(j-s  + 1)  (38) 

4 

v*"6(j)  =  '/<  E,  v*,4(j-s+i) 

i=l  ' 

R 

K*JJ(y)  =  1/8  E  V*A(j-s+ 1)  , 
which  generalizes  readily  as 

*o 

V*4n«(j)  =  1  / na  E  K^tf-s+l)  ,  (39) 

na  =  2,4,8,16 

J  -  ("o+l),  («0  +  2),...  . 

If  the  maximum  and  minimum  values  of  in-scan 
resolution  index  are  denoted,  respectively,  as  nmax  and 
nmin,  the  number  of  delay  elements  in  each  IAV  sec¬ 
tion  (Fig.  11)  is  given  by 

V  =  (vAf.)-'  •  (4°) 

For  the  example  of  Eq.  21 

A/***  *  (64/4) -1  «  15  ,  (41) 

delay  elements  per  IAV  section  are  obtained.  For  this 
example,  IAV  provides  outputs  corresponding  to  the 
following  five  in-scan  resolutions 

n  =  4,8,16,32,64  .  (42) 

The  number  of  required  delay  elements  can  become 
quite  large  for  IAV  sections  structured  as  per  Fig.  11. 
For  example,  if  instead  of  Eq.  42 

n  =  1,2,4,8,16,32,64,128,  (43) 

the  number  of  delay  elements  per  IAV  section  increases 
from  IS  to  127  in  spite  of  the  fact  that  the  number 
of  in-scan  resolutions  has  only  been  increased  from 
5  to  8. 

Depicted  in  Fig.  12  is  an  alternative  structure  for 
the  IAV  sections,  modified  from  Fig.  11  to  reduce  the 
required  number  of  delay  elements.  The  processing  ele¬ 
ment  labeled  “A”  in  Fig.  12  is  structured  like  the 
original  IAV  section,  Fig.  11.  Each  of  the  boxes  label¬ 
ed  “B”  in  Fig.  12  is  an  infinite-impulse-response  (HR) 
approximation  to  Eq.  39  for  the  appropriate  value  of 

no- 

Particularly  simple  IIR  structures,  shown  in  Figs. 
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Figure  14.  Approximate  realization  of  Hn(z),  Figure 
13.  Since  n  Is  a  power  of  2,  both  fixed  gain  amplifiers 
are  realizable  as  bit  shifts. 


Figure  12.  iAV  section,  modified  to  minimize  the 
number  of  required  delay  elements. 


exp(-1/n) 


Vk.no) 


Figure  13.  Suboptlmal  approximation  of  Hn(z)  requir¬ 
ing  just  one  memory/delay  element. 


13  and  14,  require  just  one  sample  delay  per  filter  func¬ 
tion,  Hn(z).  Following  this  approach,  the  IAV  section 
corresponding  to  the  eight  in-scan  resolutions  given 
by  Eq.  43  can  be  implemented  with  just  18  delay 
elements  rather  than  127  elements,  as  would  be  re¬ 
quired  by  the  original  approach  (Fig.  11)  alone. 

Of  course,  the  reduced  number  of  processing 
elements  achieved  via  HR  processing  (Fig.  13)  relative 
to  finite-impluse  response  (FIR)  processing  (Fig.  11) 
is  at  the  expense  of  reduced  SNR.  A  relevant  discus¬ 
sion  of  this  point  is  provided  in  [22].  The  SNR  penal- 
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Figure  15.  Alternative  IAV  structure. 


ty  involved  in  developing  IIR  approximations  to  Eq. 
39  depends  on  the  order  of  the  IIR  filter;  procedures 
for  designing  HR  structures  subject  to  design  criteria 
appropriate  to  the  present  application  are  given  in  [23, 
Chap.  8]  and  [24]. 

Yet  another  approach  to  designing  IAV  sections  to 
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minimize  the  number  of  processing  elements  is  shown 
in  Fig.  IS.  The  components  labeled  “A"  in  Fig.  IS  are 
structured  as  per  Fig.  11;  the  component  labeled  “B” 
in  Fig.  IS  is  structured  like  the  detection  preprocessor 
(Fig.  8/Eq.  23),  but  modified  for  16-fold  resolution 
reduction  and  decimation-in-time. 

The  development  of  a  multi-resolution  spatial  filter 
array  via  the  structure  of  Figs.  10  through  IS  is  an 
original  contribution  of  the  work  described  in  this 
report.  Some  virtues  of  this  approach  are  that 


•  The  image  data  are  processed  sequentially  by  col¬ 
umn;  since  scanned  linear  arrays  acquire  imagery 
sequentially  by  column,  the  processing  approach 
described  eliminates  the  need  to  buffer  large 
amounts  of  data;  and 

•  The  total  processor  comprises  a  large  number 
of  simple,  identical,  modules,  operating  in 
parallel,  and  thus  appears  inherently  well-suited 
to  real-time  applications. 


TM*  JOHN*  HOFtOM*  UNIVtMITY 

APPLIED  PHYSICS  LABORATORY 

iAUML,  MAAYIANO 


8.0  THE  PARALLEL  RANK  OF  SPATIAL  FILTERS 


The  multi-resolution  filter  bank  is  realized  by  con¬ 
catenating  XAV  and  1AV,  as  shown  in  Fig.  9.  A  single 
input /output  relationship  for  these  two  combined 
filters  is  obtained  by  substituting  Eq.  33  into  Eq.  39, 
with  the  result 


Co)  =  i,  y'^u-s+D  m 

which  may  be  written  in  a  compressed  notation  as 
S'*”°  =  V',<&rect(kln0)  ,  (45) 


where 


The  meaning  of  Eqs.  44  and  45  is  explained  with 
the  aid  of  Fig.  16.  As  shown  in  the  figure,  the  effect 
of  image  convolution  with  a  rectangular  smoothing 


1  0-n+t) 


^(l+k-1) 

Fin*  resolution  Image 
(input  image) 


Reduced 
resolution  Image 
(output  image) 


(l-k+1) 


K1-4  =  (  V)  *{J)  )u  =  image  at  output  of  detection 
preprocessor.  Fig.  8 

j/'4"0  =  (  ^0)L  =  reduced  resolution  image 
with  indices  (kt  4 na) 

reef  (A.  n0)  =  2D  rectangle  function,  of  dimen¬ 
sions  k  x  n0  and  amplitude 
( 1 /*«„)• 


Figure  16.  Convolution  with  a  rectangular  kernel  may 
be  Interpreted  ae  a  mapping  from  a  fine  resolution  Im¬ 
age  (Figure  left)  to  a  reduced  resolution  Image  (Figure 
right).  The  (l.j)  —  pixel  value  In  the  reduced  resolution 
Image  Is  obtained  as  the  average  of  the  pixel  values 
within  a  rectangular  region  (shaded  box)  of  dimensions 
k  x  n,  whose  upper-right-comer  Is  located  at  the  (IJ) 
—  pixel  location  In  the  Input  image.  As  shown  In  the 
Figure  resolution  eduction  produces  some  Image 
shrinking. 


(a)  Input  image, 
(16x32)<*ctangl* 
(M)  -  (1.1) 


(c)  A  reduced  resolution  (•)  A  reduced  resolution 

I  megs  at  output  of  IAV  Image  at  output  ol  XAV 

<k,n)  .  (1,32)  (k.n)  -  (KM) 


(g)  Reduced  resolution  image, 
convolution  kernel  matched 
to  target  dimensions 
(k.n)  *  (IS, 32) 


(b)  Sam*  data  aa  (a)  above  (d)  Sam*  data  as  (c)  above 


(f)  Sam*  dat*  as  (e)  above 


(h)  Same  data  as  (g)  above 


Figure  17.  Rectangle  test  image,  at  original  resolution,  (a)  —  (b),  and  at  3  reduced  resolutions,  (c)  —  (h).  Data 
are  prior  to  background  subtraction. 
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function  is  to  perform  a  mapping  from  a  fine  resolu¬ 
tion  image  to  a  reduced  resolution /smoothed  image, 
in  which  the  (ij)- pixel  value  in  the  reduced  resolu¬ 
tion  image  is  obtained  as  the  average  of  the  pixel  values 
within  a  rectangular  region  of  dimensions  kxn0, 
whose  upper-right-hand  corner  is  located  at  the  ( ij )- 
pixel  location  in  the  input  image. 

The  XAV/IAV  spatial  filter  array  was  implemented 
as  a  FORTRAN  program  and  exercised  against  several 
simple  test  targets  to  confirm  correct  operation.  Figure 
17  depicts  the  results  of  one  such  simulation,  show¬ 
ing  the  appearance  of  a  16x32  rectangular  test  image 
at  three  reduced  resolutions.  Detailed  numerical 
evaluations  of  similar  output  images  for  bar  and  point 
targets  (for  which  the  convolution  products  in  Eq.  45 
are  exactly  calculable)  were  used  to  verify  correct  func¬ 
tioning  of  the  FORTRAN  simulation. 

To  illustrate  the  SNR-enhancing  effect  of  spatial 
filtering,  white  gaussian  noise  was  added  to  the  rec¬ 
tangle  test  image  to  reduce  the  input  image  SNR  (target 
contrast /RMS  noise)  to  unity,  i.e., 

SNRC  =  1  . 

The  result  is  shown  as  the  upper-left  image  in  Fig.  18. 
Numerical  analysis  of  the  images  output  from  IAV  in¬ 
dicates  that  SNR  is  enhanced  predictably  as  a  func¬ 
tion  of  spatial  resolution.  In  particular,  the  (16,32)- 


filter  output,  shown  bottom-left  in  Fig.  18,  is  found 
to  have  SNR  very  nearly  equal  to  (16x32)w  = 
I6\f2.  Also  shown  for  illustrative  purposes  in  Fig.  18 
is  the  effect  of  applying  spatial  filtering  to  a  frigate 
silhouette. 


yk.n 


16x32  Rectangle  Frigate  silhouette 


Figure  18.  Multiresolution  processing  enhances 
SNR.  Four  Images  (left)  are  for  rectangle  input  Image; 
four  images  (right)  are  for  frigate  silhouette.  Images 
at  original  resolution,  having  SNR  ■  1,  are  in  the  top 
row  (rectangle  is  top-left;  frigate  is  top-right). 


34 


THE  JOHNS  hOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


9.0  BACKGROUND  ESTIMATION  AND  SUBTRACTION  (BES) 


The  need  for  background  estimation  and  subtrac¬ 
tion  (BES)  is  explained  with  the  aid  of  Fig.  19.  The 
detection  problem  is  conceptualized  as  being  one¬ 
dimensional  in  Fig.  19  to  simplify  the  discussion. 

Figure  19  depicts  a  target  of  size  yi  pixels  observ¬ 
ed  against  three  types  of  backgrounds,  viz.,  a  cons¬ 
tant  background  of  unknown  intensity  (part  a);  a 
background  of  constant  but  unknown  intensity  gra¬ 
dient  (part  b);  and  a  “slowly  varying”  background 
(part  c).  Successful  detection  requires  that  the  pro¬ 
cessor  establish  an  intensity  threshold,  T(J),  that  rides 
above  the  unknown  and  possibly  varying  background, 
B(j),  by  a  fixed  intensity  offset,  0. 

Figure  20  depicts  a  conceptual  means  for 
establishing  the  threshold  and  accomplishing  detec¬ 
tion  as  per  the  examples  in  Fig.  19.  Estimating  under 


Figure  19.  Examples  of  1-D  detection  problems  il¬ 
lustrating  the  need  for  BES.  Part  a),  target  of  dimen¬ 
sion  7s  (pixels)  in  constant  background.  Level  ‘a’  is 
unknown.  Part  b),  target  embedded  In  constant  slope 
background.  Background  parameters  'a'  and  *b’  are 
unknown.  Part  c),  target  observed  against  slowly  vary¬ 
ing  background. 


the  null  target  hypothesis  what  the  background  would 
have  been  at  each  position  j,  £(j)>  appears  to  be  an 
essential  element  in  developing  a  threshold  T(J)  for 
targe!  detection  against  spatially  varying  backgrounds. 

Equivalent  to  Fig.  20,  the  detection  process  can  be 
conceived  as  being  accomplished  in  two  stages,  viz., 
background  estimation  and  subtraction  (BES) 

C(j)  =  V{J)—B(j)  ,  (46) 

followed  by  threshold  detection 

C(J)  <  «  (47) 

A  process  analogous  to  Eq.  46  is  separately  applied 
to  all  waveforms,  i.e., 

-  Bk-"U )  ,  (48) 

as  shown  in  Fig.  9. 

A  number  of  candidate  filters  for  developing  the 
background  estimate,  B,  is  discussed  next. 

Assuming  that  the  direction  of  scan  is  parallel  to 
the  horizon,  the  background  estimate  may  be 
developed  as 


Bk-"(j)  =  Vk:"(j-nu)  .  (49) 

In  the  event  that  the  background  level  is  found  to 
contain  appreciable  variation  over  in-scan  angles  *na, 
an  alternative  background  estimator  is 


Detected 

Intensity 

VO) 


Offset,  n 


Figure  20.  Conceptual  signal  processor  for  im¬ 
plementing  threshold  detection  as  per  Figure  19 
examples. 
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(5°) 

=  [  Vk'"n(J-Snu/A)  +  Vkl  "/1(j+  5nu/4)  I  /2. 

The  image  areas  (“windows")  used  by  the 
background  estimators,  Eqs.  49  and  50,  are  depicted 
in  Fig.  21.  The  “window  gaps”  in  Fig.  21b,  separating 
the  background  windows  from  the  target  window,  are 
intended  to  prevent  target  self-thresholding  (with  a  con¬ 
sequently  degraded  SNR)  when  the  target’s  in-scan  ex¬ 
tent,  over-spills  the  target  window;  i.e.,  when 

net  <  <  2 na  .  (51) 

The  variance  of  a  random  process  X(t)  is  denoted  as 

Var  (X)  =  E  \[X-E(X)V\  >  (52) 

where  £(  • }  indicates  expectation  over  sensor  noise 
statistics. 

For  spatially  uniform  backgrounds,  specifically  ex¬ 
cluding  the  possibility  of  a  target  in  the  scene,  we  define 

«*\  =  VarlVfU)]  ,  (53) 

where  the  presumed  lack  of  dependence  of  a2k  n  on  / 
is  discussed  below. 

Since  K*’"(y)  is  statistically  independent  of  B*'a(j) 
for  all  background  estimators  of  interest,  it  follows 
from  Eqs.  48  and  53  that 


Origin  of  coordinates 


Figure  21.  Image  areas  used  by  two  background 
estimators.  Part  a)  corresponds  to  Eq.  (49);  part  b)  cor¬ 
responds  to  Eq.  (50).  Cross-channel  (x)  dimensions  are 
in  channel  units;  In-scan  (s)  dimensions  are  In  dwells. 
Scan  direction  (s)  is  parallel  to  the  horizon. 


Var[Cfn{j)]  =  a\„  +  Var[Bf'(j))  .  (54) 

It  follows  from  Eqs.  49  and  50  that,  for  both  of  these 
background  estimators, 

Var[B?(j))  =  4„  ■  (55) 

Thus,  from  Eqs.  54  and  55, 

Var[C*nU)]  =  2 a\n  ,  (56) 

i.e.,  for  both  background  estimators  considered  thus 
far,  the  differencing  process  used  to  develop  estimated 
contrast  degrades  two-fold  the  variance  of  the 
waveforms. ' 

A  relatively  less  noisy  contrast  can  be  developed  by 
using  a  largei  background  window;  e.g., 

B*”(j)  =  K*‘2"(y  — 3«m/2)  ,  (57) 

depicted  in  Fig.  22.  It  follows  from  Eq.  57  that 

Var[B*”U)}  =  4/2  •  (58) 

and  from  Eqs.  58  and  54  that 

VarlCfU))  =  1-5 o\n  .  (59) 

Comparing  Eqs.  56  and  59,  we  see  that  the  background 
estimate  provided  by  Eq.  57  results  in  a  lower-noise 
contrast  estimate. 

The  use  of  larger  background  windows  (e.g.,  Fig. 
22)  provides  improved  acquisition  performance  against 
unstructured  backgrounds  at  the  expense  of  reduced 
performance  against  structured  (“cluttered”) 
backgrounds.  The  trade-offs  involved  in  selecting  the 
MRSI  background  estimator  window  dimensions  are 
in  many  respects  analogous  to  those  encountered  in 
designing  signal  processors  for  detecting  spatially 
unresolved  targets,  commonly  known  as  infrared 
search  and  track  (IRST)  devices  [25]. 


Figure  22.  Image  area  used  by  low-noise  background 
estimator,  Eq.  (57).  Scan  direction  (s)  is  parallel  to  the 
horizon. 
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Developing  the  background  estimate  for  paralld- 
to-horizon  scanning  via  Eq.  49,  50,  or  57  imposes  an 
additional  memory  requirement  beyond  that  already 
required  by  IAV. 

Assuming  that  the  direction  of  scan  is  perpendicular 
to  the  horizon,  the  background  estimator  analogous 
to  Eq.  49  is 

B*"<j)  =  y-;k  u)  .  (60) 

Image  areas  relevant  to  Eq.  60  are  depicted  in  Fig.  23. 

Analogous  to  Eqs.  50  and  57,  one  now  has  for  cross¬ 
horizon  scanning 

B?U)  =  [  Vk'!k"(j)  +  C20')]/ 2  ,  (61) 
and 

B*nU)  =  y)k;"k  U)  ■  (62) 

Developing  background  estimates  for  cross-horizon 
scanning  via  Eq.  60,  61,  or  62  imposes  no  additional 
memory  requirement  beyond  that  already  required  by 
IAV. 

It  should  be  noted  that  the  need  for  implementing 
background  suppression  processing,  and  the  likely  ac¬ 
curacy  of  any  particular  background  estimator,  is  en¬ 
tirely  dependent  on  the  spatial  structure  of  the  ocean’s 
thermal  emission.  The  sea  radiance  is  generally  a  func¬ 
tion  of  many  variables,  including  wind  speed,  sea 
height,  nadir  angle  of  the  viewing  line-of-sight,  viewer’s 
height  above  water,  cloud  cover,  sunglints,  and  sur¬ 
face  slicks  [26]. 

In  the  processor  simulation  (App.  A),  the 
background  estimator  described  by  Eq.  50  and  Fig. 
21b,  performing  linear  trend  removal  from  the 
background,  was  chosen.  Linear  (or  nearly-linear) 
trends  may  be  caused  by  instrumental  effects  as  well 
as  by  true  thermal  gradients  in  the  scene. 

Figure  24  provides  an  example  of  a  measured  ther¬ 
mal  image  in  which  the  sea  background  appears  to 
display  a  nearly  linear  trend  with  azimuth.  It  cannot 
be  said  with  certainty  whether  the  trend  observed  in 
this  data  is  due  to  an  instrumental  effect  or  scene  ther¬ 
mal  variation.  Also,  the  magnitude  of  the  apparent 
change  in  background  level  from  one  end  of  the  scan 
to  the  other  is,  in  this  case,  still  much  smaller  than 
the  ship’s  contrast  (cf.  Fig.  32,  top). 


Continuing  the  example  that  led  to  Fig.  17,  Fig.  25 
depicts  a  rectangular  test  target  as  seen  at  four  resolu¬ 
tions,  after  background  subtraction.  The  images  in  Fig. 
25,  like  those  in  Fig.  17,  are  noise  free. 


Origin  of  coordinates 


Figure  23.  Image  area  used  by  background  estimator, 
Eq.  (60),  for  scan  direction  (s)  perpendicular  to  the 
horizon. 


100  200 
Pixel 


Figure  24.  Measured  thermal  Image  of  a  ship  observ¬ 
ed  against  a  sea  background,  In  which  a  scan  line 
through  the  background  shows  a  nearly  linear  trend. 


37 


THE  JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


(a)  (K,n)  =  (1,1)  (c)  (k,n)  =  (1,32)  (e)  (k,n)  =  (16,1)  «3)  (k,n)  =  (16,32) 


Pigurt  25.  Continuation  of  example  from  Figure  17,  depicting  rectangular  test  target  at  4  resolutions,  after 
background  subtraction. 
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10.0  THRESHOLD/DECISION  (TD)  AND  NOISE  ESTIMATION  (NE) 


Threshold/decision  processing  (TD)  currently  im¬ 
plemented  in  the  simulation  software  (App.  A)  is  based 
on  performing  a  forced-choice  detection  in  each  frame. 
Thus,  rather  than  perform  the  threshold  comparison 
indicated  in  Eq.  47,  our  TD  processor  takes  the  form 
shown  in  Fig.  26.  Each  contrast  image 

1C*'"  0)l„ 

is  normalized  with  respect  to  a  noise  standard  devia¬ 
tion,  ak  n  to  obtain  a  “signal-to-noise  image”.  The  en¬ 
tire  set  of  normalized  contrast  images  is  then  examin¬ 
ed  for  a  global  maximum  to  provide  the  forced-choice 
detection.  No  image  storage  is  required  to  implement 
TD,  since  the  search  for  the  “brightest  pixel”  is  per¬ 
formed  one  column  at  a  time.  As  a  practical  matter 

max  [C*"0)1  =  max  W<J)  ,  (63) 

ij.k.n  j 

is  calculated,  where 

W(J)  =  max  [C*"0)1  •  (64) 

i,k,n 


TD  can  be  modified  simply  to  accommodate  multi¬ 
ple  target  detections  within  a  large  field-of-view  by  im¬ 
plementing  Eq.  63  as  a  search  for  local  maxima  in 
W(j)  -  so  long  as  the  various  targets  are  spatially  well- 
separated  in  azimuth. 

The  action  of  the  noise  estimation  (NE)  component 
of  Figs.  9  and  26  is  discussed  next.  The  purpose  of 
NE  is  to  establish  estimates  for  the  noise  variances 
a\  „  defined  in  Eq.  53. 

The  root-mean-square  (RMS)  noise  levels  of  the  in¬ 
dividual  cross-scan  channels  are  denoted  as 

a,,  <  =  1,2,...,  /  , 

where  u,  is  the  RMS  value  of  ar.  additive  noise  com¬ 
ponent  of  waveform  VX(J)%  a2  corresponds  similarly 
to  y2U),  etc.  Assuming  that  the  noise  level  is  about 
the  same  for  each  channel 

Of  =  a 0  ,  /=1,2,...,/  ,  (65) 

and  that  the  additive  noise  is  signal-independent, 


The  “high  data  rate”  part  of  the  search  process  is  con¬ 
tained  in  Eq.  64.  For  the  set  of  13  spatial  resolutions 
given  by  the  example  of  Fig.  6,  and  assuming  that  the 
seeker  has  1  =  128  detector  channels  (as  per  Eq,  36), 
it  can  be  shown  that  Eq.  64  implies  that,  at  every  fourth 
sampling  instant,  TD  must  find  the  largest  of  Nw 
values,  where 

Nw  =  4/(2)  +  5/(4)  +  4/(8)  , 

with  I(k)  given  by  Eq.  34.  Thus,  implementation  of 
Eq.  64  requires  that  on  every  fourth  sampling  instant, 
TD  find  the  maximum  of  1617  data  values. 


Figure  26.  Threshold/decision  (TD)  logic  for  forced- 
choice  detection. 


da/d  F(=  0  ,  (66) 

the  RMS  noise  of  waveform  Vfn(j)  at  the  output  of 
1AV  is  given  by 

=  o0/(knV‘  .  (67) 

If  the  different  cross-scan  channels  have  appreciably 
different  noise  levels,  invalidating  Eq.  65,  it  may  then 
be  desirable  to  include  a  noise  calibration  mode  in  the 
signal  processor.  This  entails  illuminating  the  focal 
plane  uniformly  (e.g.,  by  staring  into  the  radiant  source 
used  also  for  dc-restoration),  and  then  using  the  sample 
variance  of  as  an  estimate  for  a2kn(i).  The 

RMS  noise,  akll,  in  Eq.  67  is  then  replaced  by  akn(i). 

The  possibility  of  implementing  a  noise  calibration 
mode  (as  an  alternative  to  implementing  a  “hardwired” 
calibration,  Eq.  67),  is  indicated  in  Fig.  9  as  a  dashed 
connection  between  IAV  and  NE.  Our  software  simula¬ 
tion  (App.  A)  does  not  in  fact  use  Eq.  67.  At  some 
expense  in  execution  time,  the  simulation  prior  to  each 
forced-choice  detection  executes  a  calibration  run 
against  a  noise-only  frame,  calculating  values  of  akn 
as  sample  variances,  rather  than  by  using  the  closed 
form  expression,  Eq.  67,  that  is  valid  only  for  white 
noise. 
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11.0  PROCESSOR  SIGNALTO-NOISE  MODEL 


In  this  section,  we  complete  the  simple  processor 
signal-to-noise  model  (cf.  Eqs.  1,  13,  and  19) 

r  _  *INpix  (70) 

i,A»p/R(a(3)»  , 

by  showing  that  the  choice  in  Section  5.0  of  rectangular 
response  functions  provides  a  simple  geometrical  basis 
on  which  to  calculate  the  “processor  efficiency’’,  tj. 

The  geometry  of  the  calculation  is  depicted  in  Fig. 
27,  in  terms  of  two  locally  defined  quantities,  A  and 
B,  defined  as  follows: 

A  -  number  of  target  pixels  within  a  rectangular 
region  of  integration 

B  =  number  of  background/non-target  pixels 
within  the  rectangular  region  of 
integration. 

The  contributions  of  quantities  A  and  B  relevant 
to  calculating  PSI  and  ISI  SNR  performance  are  given 
in  Thble  2.  For  example,  the  increase  in  RMS  noise  is 
always  equal  to  the  square  root  of  the  number  of  pix¬ 
els  within  the  area  of  integration,  which  for  PSI  (i.e., 
perfect)  processing  is  simply  the  ship’s  area  (defined 
in  Eq.  70  as  NP1X),  and  which  for  ISI  (i.e.,  imperfect) 
processing  is  the  area  of  the  rectangular  response 
function, 

kn  =  (A+B)  . 

The  third  column  in  Thble  2,  SNR  gain  (G),  is  ob¬ 
tained  simply  as  the  ratio  of  the  entries  in  the  first  and 
second  columns. 

The  processor  efficiency,  ij,  is  calculated  by  ratio- 
ing  the  third  column  entries  from  Table  2, 


Figure  27.  Geometry  for  calculating  processor  effi¬ 
ciency,  i),  entering  Into  Eq.  (70). 


Table  2 

Contribution  of  quantities  A  and  B,  from  Fig.  27,  relevant 
to  calculating  PSI  and  ISI  SNR  performance 


Processor 

Signal 

Enhance¬ 

ment 

RMS 

Noise 

Increase 

SNR 
Gain  (G) 

Perfect  Synchronous 
Integrator  (PSI) 
Imperfect  Synchronous 
Integrator  (ISI) 

Npix 

A 

iynx 

(kn)* 

iyPIX 

A/(kn)* 

i?  =  G (ISI) /G (PSI)  =  A/(kn  •  NPIX)‘A  .(71) 

Equation  71  is  the  desired  expression  for  processor 
efficiency,  t>,  which,  taken  together  with  Eq.  70,  com¬ 
prises  what  we  call  the  processor  SNR  model. 

We  illustrate  the  application  of  Thble  2  and  Eq.  71 
with  an  example,  in  which  it  is  assumed  that  the  target 
dimensions  are  one  pixel  in  elevation  and  m  pixels  in 
azimuth  and  that  each  of  the  m  target  pixels  has  the 
same  brightness.  Also,  we  uniquely  define  an  integer, 
n0,  in  terms  of  m,  as 

no  <  m  <  2na  ,  (72) 

where  n0  is  an  integer  power  of  2.  It  follows  from  Eq. 
72  that  the  MRSI  processor  will  estimate  the  target’s 
length  as  being  either  rt0  or  2n0  pixels,  depending  on 
which  of  the  (k,n)  =  (l,/i0)  or  (k,n)  =  (l,2n0) 
filter  output  amplitudes  is  largest.  The  entries  in  Thble 
3  are  calculated  from  the  second  row  of  Thble  2  (ISI), 
by  setting  k=  1  and  A=na  (when  n  =  n0)  or  A  =m 
(when  r,  =  2no).  The  processor  gain,  G,  is  obtained 
from  Thble  3  as 

G(m,)  =  maxG(n)  =  max[n'A0,  m/(2n0)'/,)t (73) 

n  n 

where  m,  is  an  initial  processor-derived  estimate  for 
target  size. 

It  follows  also  from  Thble  2  (first  row,  last  column) 
that 

G(PSI)  =  m'/x  .  (74) 

Thking  the  ratio  of  Eqs.  73  and  74,  we  obtain  the  pro¬ 
cessor  efficiency 
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Table  3 

Quantities  relevant  to  calculating  MRSI  processor  SNR 
gain  (G)  and  efficiency  fa)  for  the  example  discussed 
in  connection  with  Eq.  72.  Quantity  A  Is  defined  as  In 
Fig.  27. 


n 

A 

G(n)=A/n* 

»;(n)=A/(nm)'/‘ 

no 

no 

(n0/m)l/; 

2n0 

m 

m/(2n0)/1 

(m/2n0)"‘ 

V  *  » i(m,)  =  G(mt)/m'A  (75) 

where  the  quantities  ij(n)  are  given  in  Thble  3. 

This  example,  Eqs.  72  through  75,  is  developed 
numerically  in  Ikble  4  for  target  dimensions 

l  s  m  s  32, 


some  results  of  which  are  plotted  in  Figs.  28  and  29. 

The  SNR  gain,  Fig.  28,  has  a  step-stair  variation 
with  target  size,  intervals  of  constant  gain  alternating 
with  intervals  of  constant  slope.  The  breakpoints  in 
Fig.  28  may  be  calculated  by  setting  equal  the  quan¬ 
tities  G(n0)  and  G(2n0)  from  Thble  3 

nv‘0  =  m/(2n0)'A  j  ,  n0  =  1,2,4,8,16 . 

m~mr 

Thus,  the  Fig.  28  breakpoints  occur  at  non-integer 
values  of  m,  denoted  by  m,,  and  given  by 

mr  =  2"*,  r  =  1,2, 3, 4 .  (76) 

Using  vertical  bar  targets  in  Gaussian  noise,  we  have 
validated  Eq.  73  using  the  MRSI  simulation  provid¬ 
ed  in  App.  A.  In  performing  this  simulation,  we  used 
an  input  SNR  of  unity,  SNRC  -  1,  and  an  elevation 
field-of-view  of  I =64  pixels.  Estimating  the  SNR  gain, 
G,  as  the  average  value  of  64  values  obtained  for  the 
array  of  64  elevation  channels,  the  results  appear  ex¬ 
actly  as  shown  in  Fig.  28;  i.e.,  the  gain  values  obtain¬ 
ed  from  Eq.  73  and  those  obtained  via  simulation  are 
indistinguishable  on  the  scale  of  Fig.  28. 


Flgura  29.  SNR  efficiency,  from  Table  4.  Theoretical¬ 
ly  perfect  performance  corresponds  to  >>  =  1. 


Figure  28.  SNR  gain  from  Table  4.  Solid  curve,  G(PSI) 
=  m*;  dotted  values,  G(ISI)  from  Table  4.  Values  of 
breakpoints  m,  given  by  Eq.  76. 
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Table  4 

Calculation  of  MRSI  processor  gain  (G)  and  efficiency  (i?)  as  a  function  of 
target  size  (m)  for  the  example  discussed  in  connection  with  Eqs.  72  through 
75 


m 

no 

2n0 

max{  n'0'  , 

m/(2n0)«} 

=  G(ISl) 

T)  =  G/m‘; 

m. 

1 

1 

2 

1.00 

0.71 

1.00 

1.00 

1 

2 

2 

4 

1.41 

1.00 

1.41 

1.00 

2 

3 

2 

4 

1.41 

1.50 

1.50 

0.87 

4 

4 

4 

8 

2.00 

1.41 

2.00 

1.00 

4 

5 

4 

8 

2.00 

1.77 

2.00 

0.88 

4 

6 

4 

8 

2.00 

2.12 

2.12 

0.87 

8 

7 

4 

8 

2.00 

2.47 

2.47 

0.94 

8 

8 

8 

16 

2.83 

2.00 

2.83 

1.00 

8 

9 

8 

16 

2.83 

2.25 

2.83 

0.94 

8 

10 

8 

16 

2.83 

2.50 

2.83 

0.89 

8 

11 

8 

16 

2.83 

2.75 

2.83 

0.85 

8 

12 

8 

16 

2.83 

3.00 

3.00 

0.87 

16 

13 

8 

16 

2.83 

3.25 

3.25 

0.90 

16 

14 

8 

16 

2.83 

3.50 

3.50 

0.94 

16 

15 

8 

16 

2.83 

3.75 

3.75 

0.97 

16 

16 

16 

32 

4.00 

2.83 

4.00 

1.00 

16 

17 

16 

32 

4.00 

3.01 

4.00 

0.97 

16 

18 

16 

32 

4.00 

3.18 

4.00 

0.94 

16 

19 

16 

32 

4.00 

3.36 

4.00 

0.92 

16 

20 

16 

32 

4.00 

3.54 

4.00 

0.89 

16 

21 

16 

32 

4.00 

3.71 

4.00 

0.87 

16 

22 

16 

32 

4.00 

3.89 

4.00 

0.85 

16 

23 

16 

32 

4.00 

4.07 

4.07 

0.85 

32 

24 

16 

32 

4.00 

4.24 

4.24 

0.87 

32 

25 

16 

32 

4.00 

4.42 

4.42 

0.88 

32 

26 

16 

32 

4.00 

4.50 

4.50 

0.90 

32 

27 

16 

32 

4.00 

4.77 

4.77 

0.92 

32 

28 

16 

32 

4.00 

4.95 

4.95 

0.94 

32 

29 

16 

32 

4.00 

5.13 

5.13 

0.95 

32 

30 

16 

32 

4.00 

5.30 

5.30 

0.97 

32 

31 

16 

32 

4.00 

5.48 

5.48 

0.98 

32 

32 

32 

64 

5.66 

4.00 

5.66 

1.00 

32 

43 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


12.0  TARGET  SIZE  ESTIMATION 


Our  main  objective  in  developing  MRSI  was  to  find 
a  means  for  SNR  enhancement,  and,  at  least  for  the 
example  of  Section  11.0,  it  appears  that  this  has  been 
achieved.  Inspection  of  Figs.  28  and  29  shows  that  the 
SNR  gain  achieved  by  MRSI  is  quite  close  to  that 
achieved  by  the  theoretically  optimum  perfect  syn¬ 
chronous  integrator  (PSI). 

A  second  objective  was  to  have  MRSI  develop  an 
estimate  for  target  size,  and  in  this  interest,  the  cur¬ 
rent  algorithm  clearly  needs  improvement.  Inspection 
of  Figs.  30  and  31  shows  that  the  target  size  estimates 
developed  by  MRSI  can  easily  be  30%  in  error  for  the 
previously  discussed  example.  However,  as  discussed 
in  this  section,  it  appears  straightforward  to  develop 
simple  variants  of  the  original  algorithm  that  are 
capable  of  providing  greatly  improved  estimates  for 
target  size.  For  example,  we  propose  the  following 
refined  target  size  estimator,  appropriate  for  the 
preceding  one-dimensional  example: 

m  =  Int  [m]  ,  (77) 


Figure  30.  Initial  estimate  of  target  size,  m,,  from 
Table  4.  Dashed  line,  m  ■  m,  denotes  theoretically 
ideal  performance,  as  achieved  by  the  refined  target 
size  estimator,  Equations  (77)-(79),  In  the  limit  of  high 
SNR. 


where 

iog1  in  =  W  log1  m,  +  ( 1  -  W)  log1  m1  ,  (78) 


and 


W  =  -0.71  +  1.21  [G(m1)/C(m1)l  .  (79) 

Quantity  m,  in  Eqs.  78  and  79  was  previously 
defined  by  Eq.  73.  Similarly,  quantity  m2  is  defined 
implicitly  by 

G(m2)  =  min  C(n)  =  min  [n*,  m/(2rt*))  <80) 


Figure  31.  Size  estimate  percent  error,  derived  from 
mi  values  in  Table  4,  for  the  Initial  target  size  estimator. 
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The  function  Int  [•  ]  in  Eq.  77  denotes  rounding  of 
its  argument  to  the  nearest  integer  value. 

The  form  of  Eq.  78  was  suggested  by  the  require¬ 
ment  that 

log2  z  log2  m  >  log  m2  .  (81) 


Continuing  in  Thble  5  the  numerical  example  begun 
in  Ihble  4,  we  find  that  the  refined  one-dimensional 
target  size  estimator  provides  in  the  limit  of  high  SNR 
a  perfect  estimate  of  target  size;  i.e., 

lim  m  =  m  .  (84) 

SMR-  » 


The  coefficients  (-0.71,  1.21)  in  Eq.  79  were  derived  Although  it  would  be  straightforward  to  now  explore 
by  imposing  the  dual  requirements  via  simulation  the  noise  properties  of  our  refined  target 

A  ,  size  estimator,  such  as  the  size  estimate  variance,  this 

m  =  ml  ,  when  C(m,)  =  2  G(m2)  *  (82)  has  not  yet  been  done. 

The  preceding  one-dimensional  size  estimation 
an(*  algorithm,  Eqs.  77  through  79,  can  be  readily  extended 

A  „  to  the  development  of  refined  two-dimensional  size 

IV  =  Vi  ,  when  G(/n,)  =  G(m2)  .  (83)  estimates  as  well. 


Tabla  5 

Continuation  of  the  numerical  example  begun  in  Table  4.  Refined  estimate  of 
target  size,  m,  is  developed  via  Eqs.  77  through  79.  Comparison  of  left-most  and 
right-most  columns  indicates  that  m  =  m,  i.e.,  target  size  is  correctly  estimated 
when  signal-to-noise  is  high. 


m 

A 

m. 

A 

m2 

log2mj  logjiTij 

G(m,)  G(riij) 

W 

1-W 

m 

m 

1 

1 

2 

0 

1 

1.00 

0.71 

1.00 

0.00 

1.00 

1 

2 

2 

4 

1 

2 

1.41 

1.00 

1X0 

0.00 

2.00 

2 

3 

4 

2 

2 

1 

1.50 

1.41 

0.58 

0.42 

2.99 

3 

4 

4 

8 

2 

3 

2.00 

1.41 

1.00 

0.00 

4.00 

4 

5 

4 

8 

2 

3 

2.00 

1.77 

0.66 

0.34 

5.06 

5 

6 

8 

4 

3 

2 

2.12 

2.00 

0.57 

0.43 

5.94 

6 

7 

8 

4 

3 

2 

2.47 

2.00 

0.78 

0.22 

6.87 

7 

8 

8 

16 

3 

4 

2.83 

2.00 

1.00 

0.00 

8.00 

8 

9 

8 

16 

3 

4 

2.83 

2.25 

0.81 

0.19 

9.13 

9 

10 

8 

16 

3 

4 

2.83 

2.50 

0.66 

0.34 

10.13 

10 

11 

8 

16 

3 

4 

2.83 

2.75 

0.54 

0.46 

11.00 

11 

12 

16 

8 

4 

3 

3.00 

2.83 

0.57 

0.43 

11.88 

12 

13 

16 

8 

4 

3 

3.25 

2.83 

0.68 

0.32 

12.82 

13 

14 

16 

8 

4 

3 

3.50 

2.83 

0.79 

0.21 

13.83 

14 

15 

16 

8 

4 

3 

3.75 

2.83 

0.89 

0.11 

14.83 

15 

16 

16 

32 

4 

5 

4.00 

2.83 

1.00 

0.00 

16.00 

16 

17 

16 

32 

4 

5 

4.00 

3.01 

0.90 

0.10 

17.15 

17 

18 

16 

32 

4 

5 

4.00 

3.18 

0.81 

0.19 

18.25 

18 
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13.0  SHIP  IMAGE  DATA  BASES 


This  section  briefly  describes  the  three  ship  image 
data  bases  that  have  been  used  in  numerical  ex¬ 
periments  directed  toward  validating  both  the  MRSI 
processor  concept  and  a  software  implementation 
(App.  A).  The  three  sources  of  ship  imagery  used  thus 
far  include: 

•  NWC  FLIR  data, 

•  Mini  FLIR  data,  and 

•  APL  model  silhouettes. 

The  NWC  data  comprise  7537  thermal  images  of 
eight  different  ship  classes,  observed  from  the  air  at 
shallow  depression  angles,  at  aspects  of  30°  and  90° 
off-bow,  for  distances  of  from  6  to  20  nmi.  As  the  data 
were  obtained  during  a  time  of  very  good  atmospheric 
visibility,  the  images  consequently  have  very  high  SNR 
and  are  not  directly  suitable  for  realistically  testing  the 
limits  of  signal  processing  methods  for  SNR  enhance¬ 
ment.  For  these  data  to  be  used,  they  were  first  degrad¬ 
ed  with  progressively  increasing  levels  of  additive  noise, 
as  illustrated  in  Fig.  32. 

Of  course,  much  of  the  interest  in  working  with  these 
data,  stemming  from  the  “real  world  flavor”  of 
measured  imagery,  is  lost  when  the  data  are  syn- 


Figure  32.  High  SNR  Image  (top)  measured  under 
unusually  good  atmosphere  conditions  Is  degraded 
with  progressively  Increasing  levels  of  additive  noise 
to  simulate  observations  obtained  during  less 
favorable  weather.  To  the  right  of  each  Image  Is  shown 
a  scan  line  obtained  at  the  elevation  Indicated  by  a 
horizontal  arrow. 


Figure  33.  Measured  ship  thermal  Image  threshold- 
ed  to  display  what  appear  to  be  ac-coupllng  artifacts 
In  front  of  and  behind  the  ship,  near  the  waterline 


thetically  degraded.  Since  we  had  very  little  informa¬ 
tion  about  either  the  sensor  or  the  data  reduction  pro¬ 
cess  at  the  time  this  work  was  performed,  and  with 
what  appears  to  be  AC-ccupling  artifacts  in  some  of 
the  data  (e.g.,  Fig.  33),  we  have  made  only  limited  use 
of  these  data  thus  far  (cf.  Section  14.1).  It  is  our 
understanding  that  a  second  NWC  ship  image  data 
set  has  been  developed,  containing  images  measured 
at  lower  SNR.  Unfortunately,  the  low  SNR  NWC  data 
were  not  available  at  the  time  this  work  was  performed. 


Figure  34a.  Ship  IR  Image  digitized  from  minlFUR 
analog  video. 


Figure  34b.  Scan  line  through  the  peak  intensity  pix¬ 

el  In  the  mlnIFLIR  image  shown  in  Part  a. 
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A  second  source  of  ship  thermal  imagery  was  some 
mini  FLIR  analog  video  that  became  available  via  a 
previous  APL  program.  Several  frames  of  this  data 
were  digitized  by  K.  Constantikes  (APL/FIF)  for  our 
use  Since  a  number  of  important  sensor  characteristics 
and  the  environmental  conditions  prevailing  at  the  time 
of  measurement  were  unfortunately  not  known,  we 
have  made  only  limited  use  of  the  mini  FLIR  data  (cf. 
Section  14.1).  Visual  inspection  of  the  data,  both 
analog  and  digital,  shows  the  SNR  to  be  only  marginal¬ 
ly  adequate  for  visual  detection  (cf.  Fig.  34).  These 
data  are  therefore  of  at  least  qualitative  interest,  since, 
once  digitized,  they  can  be  used  to  exercise  low  SNR 
detection  processing  without  first  requiring  synthetic 
noise  insertion. 

Our  third  source  of  ship  image  data,  the  APL  model 
silhouette  data  base,  was  originally  developed  in  con¬ 
nection  with  an  earlier  APL  Independent  Research  and 
Development  (IRAD)  effort  [3,28].  The  data  were  ac¬ 
quired  by  digitizing  TV  images  of  scale  models  of  five 
ship  types  at  21  distinct  aspects  around  the  starboard 
bow  quadrant  for  each  of  three  simulated  ranges  (cf. 
Fig.  33).  Thus,  a  total  of  313  silhouette  images  is 
available  (3  ship  types  x  3  ranges  x  21  aspects).  The 
models  were  cast  in  lead  at  a  1:1230  scale,  and  are 
representative  of  ships  at  the  waterline.  The  images 
are  noise-free,  black/white  silhouettes,  digitized  to  312 
x  312  pixels  per  image. 

A  wide  variety  of  effects,  some  or  all  of  which  may 
figure  importantly  in  measured  imagery,  are  obvious¬ 
ly  not  present  in  the  model  silhouette  data,  viz., 


Figure  35.  APL  modal  silhouettes  are  available  for  5 
ship  types,  each  observed  at  21  aspect  angles  and  3 
ranges  [28]. 


•  background  structure  (e.g.,  due  to  waves, 
reflected  clouds,  natural  ocean  thermal  varia¬ 
tions,  and  slicks  [26]), 

•  ship  thermal  structure  (e.g.,  stack  hot  spots), 

•  atmospherics, 

•  ship’s  reflection  in  the  water, 

•  ship’s  wake,  and 

•  assorted  sensor  artifacts  (e.g.,  AC-coupling, 
“striping”  due  to  uncompensated  detector  non- 
uniformities,  vignetting,  etc). 

One  consequence  of  these  limitations,  for  example,  is 
that  the  model  silhouette  data  cannot  be  used  to  assess 
the  potential  usefulness  of  thermal  gray  scale  infor¬ 
mation  in  assisting  classification. 

Certain  real-world  characteristics  can  be  imparted 
to  the  model  silhouettes  via  modeling.  Fbr  example, 
Fig.  36  illustrates  how  model  silhouettes  (three  images 
at  figure  upper  left)  may  be  combined  with  ship  IR 
signature  data  (curve  at  figure  upper  right)  that  in¬ 
corporate  both  ship  thermal  modeling  [20]  and  a 
statistical  treatment  of  weather  effects  [12].  Seeker  sen¬ 
sitivity  is  described  by  the  NEAT  parameter  (curve  at 
figure  lower  right).  Output  from  the  image  model  is 
a  set  of  images  (figure  lower  left)  having  SNRs  that 
reflect  a  wide  variety  of  phenomenological 
considerations. 

It  is  pointed  out  in  [3]  that  the  model  silhouettes 
appear  well-suited  to  developing  ship  recognition 
techniques  based  on  ship  profile  information,  since 
the  ship  profile  may  be  substantially  the  same  for  1R 
imagery  as  for  TV  imagery.  Of  greater  interest  to  the 
present  application,  it  is  noted  that  the  well-controlled 
measurement  conditions  for  the  model  silhouettes 
facilitate  the  development  of  curve  fits  for  ship  pro¬ 
jected  area  Vs.  viewing  aspect.  It  is  found  that,  with 
relatively  small  error  (cf.  Fig.  37) 

,  ...  a  4750  sin  $  +  645  cos  6  ,  carrier  ... 

880  sin  6  +  90  cos  9  ,  frigate  , 

over  the  range  7°  <  $°  s  90°  (beam  aspect),  and 


A*,  <*)  “  -W70)  •  °°  *  *  *  7°  .  (86) 

It  follows  from  Eq.  85  that  the  processor  SNR 
model,  Eq.  70,  can  now  be  written  as 

G(B  R)  »  OT  2n/Rj  [(sine  +  0.14  cos0)/a8]  carrier 
(\6i)/R0)‘Hsin6  +  0.10  cos6)/aB]  '*,  frigate  , 
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0.000  K  02  K  0.6  K 


SNR  -  (AT7NEAT) 


Flgura  SB.  Image  modal  combines  Idealized  silhouettes  (3  Images,  upper  left)  and  ship  IR  signature  data  The 
IR  signature  (curve,  upper  right)  Incorporates  both  ship  thermal  modeling  and  a  statistical  treatment  of  weather 
effects  (121- 


over  the  range  7°  <  6  s  90°,  and,  at  smaller  aspects 

G(0,Ko)  at  G(7°,  Ro)  ,  0'  S  «  S  7V  (88) 

We  note  that  R0  in  Eq.  87  is  range  in  nmi  units,  and 
(a,  0)  are  in  mr  units. 


For  illustrative  purposes,  we  plot  in  Fig.  38  the  varia¬ 
tion  of  processing  gain  with  range,  as  obtained  from 
Eq.  87,  for  two  aspects  (90°  and  43  °),  against  a  frigate 
target.  In  obtaining  Fig.  38,  a  spatial  resolution  of  a 
=  0  =  0.15  mr  and  a  processor  efficiency  of  ij  *  1 
are  assumed. 
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Figure  37.  Ship  projected  area  vs.  viewing  angle  ob¬ 
tained  from  APL  model  silhouette  imagery  ot  a  frigate 
(solid  line).  Curve  fit  (dashed  line)  was  obtained 
heuristically. 


Figure  38.  Processing  gain  (G)  vs.  range,  against 
frigate  target,  from  Eq.  (87).  Assumed  spatial  resolu¬ 
tion  is  a  =  0  =  0.15  mr;  assumed  processor  efficien¬ 
cy  iS  t]  =  1. 
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14.0  RESULTS  OF  IMAGE-BASED  SIMULATIONS 


An  image-based  simulation  of  the  MRSI  processor 
has  been  coded  in  FORTRAN  (App.  A).  The  simula¬ 
tion  operates  on  images  of  64  x  2S6  pixels,  at  35  spatial 
resolutions.  The  program  runs*  with  minor  changes, 
on  both  the  APL  NAS  mainframe  and  on  an  IBM  PC 
AT  computer  and  requires  about  450  kbytes  of 
memory.  The  listing  provided  in  App.  A  is  for  the  IBM 
PC  AT  version  of  the  code.  All  simulation  results  pro¬ 
vided  in  this  section  (the  results  of  more  than  2800 
processed  images)  were  developed  on  an  IBM  PC  AT. 


14.1  INITIAL  QUALITATIVE  RESULTS 

Our  first  MRSI  simulations  against  ship  imagery 
were  performed  using  NWC  imagery,  synthetically 
degraded  in  SNR,  as  discussed  in  connection  with  Fig. 
32.  It  was  observed  for  these  first  simulations  that  cor¬ 
rect  detections  and  reasonably  accurate  size  estimates 
were  developed  with  input  SNR  values  of  0.2.  One  of 
the  results  of  these  early  efforts  is  shown  in  Fig.  39, 
in  which  a  ship  at  13.5  nmi  range  is  correctly  detected 
and  sized,  for  SNRC  =  0.2. 


Original  image, 
San  Diego  locale 


SNR  =  0.2 
Simulated  adverse 
atmosphere 


Forced -choice 
detection 


Figure  39.  Initial  qualitative  evidence  of  correct 
simulation  performance.  Original  high-SNR  NWC  ship 
image  (top)  was  degraded  to  SNR  =  0.2  (center)  to 
simulate  observation  in  a  less  favorable  atmosphere. 
Degraded  image  was  Input  to  FORTRAN  simulation 
(App.  A)  and  the  ship  subsequently  detected  and  siz¬ 
ed  (bottom).  Range  to  ship  Is  13.5  nml. 


A  similar  qualitative  result  was  obtained  by  exer¬ 
cising  MRSI  against  a  mini  FLIR  image  of  uncertain 
SNR  (cf.  Fig.  34).  As  shown  in  Fig.  40,  in  this  case, 
we  ran  MRSI  in  an  experimental  multi-pass  configura¬ 
tion,  and,  after  three  passes,  were  provided  by  MRSI 
with  a  highly  ship-like  composite.  The  ship  in  this  data 
was  observed  at  a  range  of  about  10  nmi. 


Figure  40.  Multipass  detection  of  low  contrast 
mlniFLIR  Image  provides  shape  information.  Original 
Image  shown  at  top;  composite  product  of  3-pass 
detection  shown  at  bottom.  Range  to  ship  is  10  nmi. 


14.2  PROCESSOR  EFFICIENCY 

The  input  images  for  the  simulations  described  in 
this  section  were  various  APL  model  silhouettes  (cf. 
Section  13.0),  to  which  was  added  a  small  amount  of 
Gaussian  noise  to  obtain 

SNRC  =  5  .  (89) 

As  discussed  in  connection  with  Fig.  26,  the  detected 
targets  are  characterized  by  MRSI  in  .erms  of  five 
parameters,  viz., 

•  position  (i0,  j0), 

•  size  (k0,  n0),  and 

•  contrast  C0. 
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In  addition,  as  discussed  in  connection  with  Fig.  26 
and  Eq.  67,  the  processor  develops  an  estimate  of  the 
noise  variance,  crkn,  for  each  spatial  filter  output.  An 
estimate  for  the  detection  SNR,  defined  in  Eq.  12,  is 
then  obtained  from  the  results  of  simulation  as 

.  (90) 

However,  from  Eq.  13 

V  =  G  •  N'p,x  =  ( SNRd„/SNRc )  •  N'plx  .  (91) 

We  obtain  from  Eqs.  90  and  91  the  following  estimate, 
17, ,  for  processor  efficiency 

Vi  =  ( SNRde/SNRc )  ■  N'plx  =  6  ■  N*x,  (92) 

where  SNRC  (Eq.  89)  and  NPjX  are  known  a  priori, 
and  is  an  output  of  the  simulation. 

Since  Eq.  90  provides  only  an  estimate  for  SNRdel, 
we  improve  the  estimate  by  performing  each  simulated 
detection  a  number  of  times  ('  ich  time  with  a  different 
random  number  seed)  and  then  averaging  over  the 
resulting  individual  SNRdel  estimates. 

Alternatively,  Eq.  71  and  the  position  and  size 
estimates  provided  by  the  simulation,  (ic,  jj  and  (k0, 
nD),  can  be  used  to  estimate  ij  as 

V2  =  A/(kono  ■  NPIX)»  .  (93) 

Calculations  of  processor  efficiency  have  been  per¬ 
formed  for  13  model  silhouettes  (Thble  6).  The  first 
series  of  Monte  Carlo  simulations  was  performed  with 
a  very  large  number  (120)  of  trials  to  obtain  an  estimate 
for  how  many  Monte  Carlo  trials  were  needed  to 
achieve  adequate  convergence  of  the  estimate  for 
SNRdel.  The  silhouette  image  used  in  all  120  trials 
was  a  frigate  profile,  for  range  =  10  nmi,  aspect  = 
90°;  the  only  differentiating  aspect  of  the  120  input 


Table  6 

Calculations  of  processor  efficiency  »>  have  been  per¬ 
formed  for  13  model  silhouettes. 


Ship 

Aspect  (°) 

Ranges  (nmi) 

Frigate 

10 

5,  10 

45 

5,  10,  15 

90 

10,  20 

Carrier 

10 

10,  20 

45 

20,  30 

90 

20,  30 

images  was  in  thr  random  number  seed  used  in 
generating  the  additive  noise.  The  results  of  these  120 
trials  are  given  in  Fig.  41  as  a  plot  of 

G  =  $NRdtl  /5  , 

as  a  function  of  the  number  of  Monte  Carlo  trials. 


Figure  41.  Estimated  processing  gain,  as  a  function 
of  the  number  of  Monte  Carlo  trials.  The  same 
silhouette  image  was  used  in  all  120  trials  (frigate  at 
10  nmi,  aspect  =  90*).  120  distinct  images  were 
created  by  adding  to  the  noise-free  ship  image  120  dif¬ 
ferent  noise  Images. 


We  surmised  from  this  numerical  experiment  that 
no  more  than  15  Monte  Carlo  trails  are  required  to 
develop  an  estimate  for  G  accurate  to  within  5  or  10% 
of  the  fully  converged  value:  All  consequent  determina¬ 
tions  of  t)  via  simulation  for  the  remaining  12  ship  pro¬ 
files  listed  in  Thble  6  were  established  as  averages  over 
15  Monte  Carlo  trials. 

Estimates  of  processor  efficiency  for  the  13  ship  pro¬ 
files,  obtained  using  Eq.  92,  are  provided  in  Thble  7. 
Estimates  of  7?  based  on  Eq.  93  are  given  in  Thble  8. 

The  dispersion  of  rj ,  values  in  Thble  7  is  bounded 
by 

|  i),  -  0.82  |  <  0.08  ,  (94) 

with  mean  and  variance 

£(<?,)  =  0.81  (95) 

Vartft)  =  (0.048)2  . 

The  dispersion  of  t;2  values  in  Thble  8  is  bounded  by 
|  rj2  -  0.76  |  <  0.05  ,  (96) 
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Table  7 

Processor  efficiency,  i»,  is  estimated  using  Eq.  92  and  results  from  Image-based 
simulations,  for  13  ship  silhouettes.  Each  value  of  6  in  this  table  is  established 
as  an  average  over  15  Monte  Carlo  Iterations;  thus,  a  total  of  195  images  was  pro¬ 
cessed  to  obtain  these  results.  Assumed  IFOV  =  0.073  mr 


Ship 

Aspect  (°) 

Ranges  (nmi) 

NPIX 

6 

i 

Frigate 

10 

5 

479 

16.77 

0.77 

10 

139 

9.32 

0.79 

45 

5 

1247 

30.69 

0.87 

10 

355 

15.66 

0.83 

15 

139 

9.49 

0.80 

90 

10 

477 

19.63 

0.90 

20 

122 

8.90 

0.80 

Carrier 

10 

10 

702 

20.66 

0.78 

20 

175 

10.00 

0.76 

45 

20 

434 

18.14 

0.87 

30 

182 

10.66 

0.79 

90 

20 

555 

20.24 

0.86 

30 

248 

11.82 

0.75 

Table  8 

Processor  efficiency  for  13  ship  profiles  as  estimated  from  Eq.  93 


Ship 

Aspect  (°) 

Ranges  (nmi) 

G(PSI) 

‘^pix 

G(ISI) 
A/(k0  nor 

A 

Frigate 

10 

5 

21.89 

16.62 

0.759 

10 

11.79 

9.19 

0.780 

45 

5 

35.31 

26.66 

0.755 

10 

18.84 

13.74 

0.730 

15 

11.79 

9.19 

0.780 

90 

10 

21.84 

17.50 

0.801 

20 

11.05 

8.75 

0.792 

Carrier 

10 

10 

26.05 

19.49 

0.748 

20 

13.23 

9.72 

0.735 

45 

20 

20.83 

16.18 

0.777 

30 

13.49 

10.25 

0.760 

90 

20 

23.56 

18.03 

0.765 

30 

15.75 

11.23 

0.713 

with  mean  and  variance 

E(ri2)  =  0.76  (97) 

Var{ $,)  =  (0.025  )*  . 

Based  on  a  total  of  195  Monte  Carlo  simulations 
with  13  different  ship  profiles,  the  efficiency  of  MRS1 
can  be  characterized  as 


t?  -  0.8  ,  (98) 

accurate  to  about  ±10%,  for  frigate  and  aircraft  car¬ 
rier  profiles,  observed  at  ranges  from  5  to  30  nmi,  and 
for  aspect  angles  from  10°  to  90°. 

The  rather  remarkable  invariance  of  ij  with  respect 
to  ship  type,  scale,  and  viewing  aspect  was  perhaps 
presaged  by  our  earlier  analytically  derived  results  for 
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bar  targets,  Fig.  29.  In  a  speculative  vein,  we  note  that 
the  bar  target  results  can  be  summarized  as 

|  v  ~  0.92  |  <  0.08  . 

Heuristically  extending  the  bar  target  results  to  two 
dimensions 

ij  as  (0.92)2  fit  0.85  , 
we  get  a  result  fairly  close  to  Eq.  98. 


14.3  THRESHOLD  SNR  FOR 
FORCED  CHOICE  DETECTION 

We  explored  systematically  the  variation  of  pro¬ 
bability  of  forced  choice  detection  (PD)  with  SNR  for 
three  ship  profiles,  involving  2500  Monte  Carlo  trials 
(Thble  9).  For  each  ship  profile,  for  each  SNR  value, 
100  Monte  Carlo  trials  were  performed;  the  relative 
number  of  times  the  declared  target  overlapped  the 
true  target  was  used  as  an  estimate  for  PD. 

Our  results  are  presented  in  Fig.  42  as  plots  of  PD 
vs.  SNR,  parametric  in  viewing  aspect.  Each  of  the 
25  data  points  in  Fig.  42  was  established  via  100  image- 
based  simulations. 

The  data  are  replotted  in  Fig.  43  as  PD  vs.  SNRdet, 
where  SNRdet  is  given  by  Eq.  12  as 

SNRdtl  m  r,  •  SNRC  •  N*x 

with  an  assumed  value  of  »;  =  0.81.  Also  shown  as 
a  solid  line  in  Fig.  43  is  a  plot  of  the  function 

PD(SNRdtl)  =  *(SNRdtl  -  2.8)  ,  (99) 

where 

$(t)  =  (l/vHZir)  J"^  exp(-t2/2)dt  . 

Interestingly,  it  has  been  found  in  psychovisual  ex¬ 
periments  performed  with  human  observers  of  elec¬ 
tronic  displays  that  a  50%  probability  of  detection 
against  rectangular  target  images  requires  a  “display 
SNR”  of  2.8  [4,  p.  86],  Equation  99  is,  as  per  the 
discussion  in  [3,  p.  197],  a  curve  fit  to  experimental 
data  descriptive  of  human  vision  system  (HVS) 
performance. 

The  fact  that  the  PD  data  for  all  three  aspect  angles 
are  nicely  fit  in  Fig.  43  by  a  single  universal  curve  sug¬ 
gests  that  PD  statistics  may  be  derived  from  the 
universal  curve  for  a  wide  variety  of  ship  target  pro¬ 


files.  "testing  the  generality  of  this  hypothesis  would, 
however,  clearly  require  more  simulation  and  analysis 
than  presented  in  this  report. 


Figure  42.  Probability  of  detection  (Pq)  vs.  channel 
SNR  (SNRC).  Frigate  silhouette  having  range  resolu¬ 
tion  product  s  0.73  nmi-mr.  Curves  are  parametric  in 
viewing  aspect:  •  *  90°,  ■  =  45°,  □  =  10*. 


Figure  43.  Detection  probability  (Pd)  vs.  detection 
SNR  (SNRqet,  defined  In  Eq.  12).  Frigate  silhouette, 
(range  x  resolution)  =  0.73  nmi-mr.  Three  symbol  types 
correspond  to  different  viewing  aspects:  •  =  90°, 
□  =  45°,  ■  =  10°.  Solid  line  is  an  analytical  curve 
fit  to  experimental  psychovisual  data,  Eq.  (99). 
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Table  9 

Probability  of  detection  statistics  were  accumulated  as  a  function  of  SNR  for  a  total 
of  2500  Monte  Carlo  trials.  Range  -  resolution  product  =  0.73  nmi  •  mr 


Ship 

Aspect  (°) 

NPIX 

Range  of 
SNR  Values 

SNR 

Increment 

Monte  Carlo  Ttials 
Per  SNR  Value 

Frigate 

10 

139 

0.15-0.60 

0.05 

100 

45 

355 

0.15-0.60 

0.05 

100 

90 

477 

0.20-0.60 

0.10 

100 
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15.0  SYSTEM  CONCEPTS  FOR  CUED  DETECTION  AND  ATR  J 


As  discussed  in  Section  1.0  the  current  study  was 
motivated  by  an  interest  in  the  antiship  missile  (ASM) 
application.  In  this  section,  we  discuss  some  top-level 
ASM  system  concepts  in  which  the  MRSI  processing 
approach  could  be  employed. 

Figure  44  depicts  a  tactical  scenario  in  which  an 
ASM  has  been  launched  in  the  general  direction  of 
a  hostile  surface  force,  with  the  problem  of  target  selec¬ 
tion  to  be  addressed  post-launch.  The  missile  seeker 
acquires  a  first  candidate  ship  target,  which  is  subse¬ 
quently  classified  as  ineligible  for  attack.  The  ASM 
continues  its  flight,  and,  as  depicted  in  Fig.  44,  ac¬ 
quires  a  second  ship,  which  it  classifies  and  engages 
as  a  high-priority  hostile  combatant. 

The  initial  pre-acquisition  surveillance  performed 
by  the  ASM  seeker  is  presumably  performed  over  a 
field-of-view  much  larger  than  that  of  typical  thermal 
imaging  systems  and  TVs.  The  initial  target  acquisi¬ 
tion  could  be  performed  by  MRSI  processing,  which 
then  directs  the  seeker’s  further  attention  to  the 
neighborhood  of  “objects”  requiring  classification. 
Thus,  the  MRSI  processor  is  followed  by  a  classifica¬ 
tion  processor  that  makes  a  series  of  judgments  regar¬ 
ding  the  nature  of  each  object,  viz.,  ship/non-ship; 
combatant/non-combatant;  hostile/friendly;  high 
value/low  value;  etc.  TWo  system  concepts  for 
classification  will  be  described. 

Figure  44  shows  a  system  concept  in  which  video 
is  telemetered  from  the  ASM  to  an  aircraft,  where 
target  classification  is  performed  and  attack  sequen¬ 


cing  initiated  by  a  human  operator.  According  to  this 
concept,  MRSI  plays  the  role  of  a  target  “cuer”,  reduc¬ 
ing  the  bandwidth  requirements  of  the  video  link  and 
unburdening  the  human  image  interpreter  from  the 
need  to  perform  visual  search  over  a  v.  ide  field-of-view. 

Figure  45  depicts  an  alternative  system  concept  in 
which  classification  is  performed  on-board  the  missile 
by  an  autonomous  target  recognition  (ATR)  algorithm. 

A  potential  problem  with  both  system  concepts  is 
that  the  seeker’s  acquisition  range  is  iikely  to  be  much 
greater  than  the  range  at  which  aijcurate  classification 
can  be  performed,  since  classification  both  by  eye  and 
by  ATR  is  likely  to  have  far  more  demanding  SNR  re¬ 
quirements  than  the  MRSI  acquisition  algorithm. 
Although  the  ASM  could  simply  flv  toward  each 
potential  target  until  the  SNR  grows  large  enough  to 
permit  classification,  this  approach  could  result  in  the 
waste  of  considerable  fuel  capacity. 

A  second,  more  satisfactory,  solution  to  the  acquisi- 
tion/classification  SNR  mismatch  problem  has  been 
suggested  by  W.  J.  Tropf  (APL).  Following  MRSI  ac¬ 
quisition,  the  seeker  could  be  directed  into  a  “classifica¬ 
tion  mode”,  in  which  scan  is  performed  very  slowly 
over  a  very  narrow  field-of-view  in  the  neighborhood 
of  the  potential  target.  The  reduced  angular  rate-of- 
scan,  combined  with  a  proportionally  reduced  post- 
detector-amplifier  electrical  bandwidth,  then  provides 
the  SNR  boost  required  for  classification,  without 
sacrificing  spatial  resolution. 

For  example,  the  classification  mode  may  be  per- 


On-board  processing 
Sensed  Image 
Reference  Image  1^^  _J 


ReacquIre/claasify/attacK 


Commerclal 

Acqulre/classlfy/dlsengage 


Figure  44.  ASM  attack  scenario,  depicting  system 
concept  In  which  target  acquisition  is  performed  by 
MRS),  and  target  classification  Is  performed  off-board 
by  a  person  inspecting  cued  images  transmitted  via 
video  link  from  the  ASM. 


Figure  45.  ASM  attack  scenario,  depicting  system 
concept  in  which  target  acquisition  Is  performed  by 
MRSI  and  target  classification  Is  performed  by  on¬ 
board  ATR. 


57 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


formed  with  an  angular  rate-of-scan  25  times  slower 
than  the  acquisition  mode, 

bc  =  d/25  , 

where  $  was  defined  previously  in  connection  with  Eqs. 
25  and  27.  If  the  detected  waveforms  are  sampled  at 
the  same  rate  during  slow  scan  (classification  mode) 
as  during  fast  scan  (acquisition  mode),  the  data  dur¬ 
ing  slow  scan  will  be  over-sampled  25-fold,  from  the 


standpoint  of  satisfying  Nyquist’s  criterion  and  preser¬ 
ving  spatial  detail.  Consequently,  the  A/D  process  can 
be  followed  by  a  pre-classification  processor  con¬ 
figured  as  per  Fig.  8  and  Eq.  23,  in  which  the  data 
are  smoothed  and  decimated-in-time  25-fold,  with  a 
consequent  5-fold  SNR  gain. 

The  cued  slow-scan  classification  mode  is  potentially 
useful  to  both  system  concepts,  Figs.  44  and  45,  i.e., 
regardless  of  whether  classification  is  performed  by 
ATR  or  by  a  human. 
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16.0  A  NOVEL  METHOD  FOR  CLOUD  CLUTTER  SUPPRESSION 


Cloud  reflections  in  the  water  are  the  most  intense 
natural  source  of  sea  surface  thermal  structure  [26], 
and  are  consequently  expected  to  offer  the  most  severe 
source  of  false  alarms  to  this  (and  other)  ship  1R  detec¬ 
tion  schemes.  In  this  section  we  present  a  no-  el  method 
for  suppressing  false  alarms  induced  by  cloud  reflec¬ 
tions  from  the  sea. 

Our  new  clutter  suppression  approach  entails  cor¬ 
elating  MRSI  detections  obtained  below  the  horizon 
(due  either  to  a  ship  or  a  cloud  reflection)  with  detec¬ 
tions  obtained  above  the  horizon  at  the  same  azimuth: 
the  presence  of  cloud  structure  above  the  candidate 
target’s  location  indicates  a  likely  false  target,  while 
the  absence  of  clouds  indicates  a  likely  true  target. 

We  illustrate  the  method  with  the  aid  of  Fig.  46, 
depicting  the  spatial  variation  of  sky  radiation 
temperature  (top)  and  sea  surface  radiometric 
temperature  (bottom),  as  measured  at  10.6  /un,  in  the 
presence  of  clouds.  These  data  were  measured  by 
JHU/APL  with  a  calibrated  dual-band  1R  radiometer 
installed  in  a  Navy  P-3  aircraft  [26].  For  example,  the 
sea  surface  thermal  artifact  labelled  A'  in  the  lower 
trace  is  clearly  indicated  as  a  cloud  reflection,  since 
it  correlates  directly  with  thermal  structure  observed 
above  the  horizon,  i.e.,  the  artifact  labelled  A  in  the 
upper  trace  of  the  figure. 

The  vector  of  attributes  attached  to  each  MRSI 
detection  (cf.  Fig.  26)  can  be  augmented  with  a  “quali- 


Flgure  46.  Spatial  variation  of  sky  radiometric 
temperature  (top)  and  sea  surface  temperature  (bottom) 
measured  by  APL  at  10.6  pm,  In  the  presence  of  clouds 
[26].  Sea  surface  artifacts  labelled  with  primed  letters 
In  the  lower  trace  are  reflections  of  cloud  structures 
labelled  with  corresponding  unprlmed  letters  In  the  up¬ 
per  trace. 


ty”  measure  derived  from  correlation  with  the  sky 
measurements. 

The  mechanization  of  our  cloud  clutter  discrimi¬ 
nant  is  simplified  by  making  use  of  the  fact  that  a 
cloud’s  contrast  against  the  sky  is  generally  several  hun¬ 
dred  times  greater  than  the  contrast  of  the  cloud’s 
reflection  against  the  sea.  For  example,  we  see  from 
Fig.  46  that  cloud/sky  contrast  is  about  50  K,  while 
the  reflection/sea  contrast  is  only  about  0.1  K.  In  ad¬ 
dition,  since  there  is  no  particular  benefit  to  be  deriv¬ 
ed  from  spatially  resolving  the  fine  details  within  the 
cloud  structure,  the  sky  radiance  data  may  be  acquired 
by  means  having  greatly  reduced  resolution  relative  to 
the  means  used  for  detecting  ships  against  the  ocean 
background.  In  Table  10  we  compare  approximate  sen¬ 
sitivity  and  spatial  resolution  requirements  for  cloud 
detection  with  corresponding  requirements  for  ship 
detection. 

As  a  first  implementation  possibility,  we  consider 
obtaining  sea/ship  radiance  data  on  a  forward  scan 
in  azimuth,  and  sky/cloud  radiance  data  on  the 
backscan.  In  this  case,  the  backscan  is  stepped  upwards 
in  elevation,  and  executed  at  a  greatly  increased  angular 
rate-of-scan,  relative  to  the  forward  scan.  Considerable 
image  smear  is  tolerable  on  the  backscan,  due  to  the 
relaxed  spatial  resolution  requirements  for  the  sky 
measurements  (cf.  Table  10).  Moreover,  if  necessary, 
excessive  smear  can  be  prevented  simply  by  widening 
the  post-detector  electronic  passband;  the  consequent 
increase  in  noise  is  accomodated  by  the  relaxed  sen¬ 
sitivity  requirements  during  the  backscan  (cf.  Ihble  10). 
Relaxed  spatial  resolution  requirements  in  elevation 
can  be  exploited  as  a  reduction  in  signal  processing 
requirements  during  the  backscan.  This  is  accomplish¬ 
ed  by  adding  together  the  detector  outputs  in  groups 
(say,  of  five  each),  to  obtain  a  single  reduced-resolution 
video  channel  for  each  group. 


Table  10 

Sensitivity  and  spatial  resolution  requirements  for 
ship  detection  and  cloud/clutter  detection  (approx¬ 
imate  values). 


Ship  Detection 

Cloud  Detection 

Field  of  Search 

Below  horizon 

Above  horizon 

Sensitivity,  NEAT 

0.1  K 

5  K 

Spatial  Resolution 

0.2  mr,  square 

1  mr,  square 
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A  second  implementation  possibility  is  to  perform 
the  sea  and  sky  measurements  simultaneously  (i.e.,  on 
the  same  scan)  by  having  a  partitioned  focal  plane: 
small  elements  at  the  bottom  of  the  focal  plane  are 
used  for  below-horizon  search,  and  much  larger  detec¬ 
tors  at  the  top  of  the  focal  plane  are  used  for  above¬ 
horizon  search.  The  optics  design  for  this  sensor  con¬ 
cept  is  complicated  by  the  requirement  for  wide  eleva¬ 
tion  instantaneous  field-of-view.  However,  this  problem 
is  partially  ameliorated  by  the  fact  that  a  considerable 
degradation  in  the  optics  blur  can  be  tolerated  over 


the  upper  part  of  the  focal  plane. 

Yet  a  third  implementation  possibility  is  to  perform 
the  sky  search  with  a  separate  telescope  The  telescope 
used  for  sky  search  can  be  very  small,  perhaps  just 
0.5-1  inch  in  aperture. 

Finally,  additional  sensor  concepts  for  nearly 
simultaneous  sky/sea  search  can  be  based  on  cross- 
horizon  scanning  (discussed  briefly  in  connection  with 
Eqs.  60-62).  The  principal  penalty  imposed  by  this  ap¬ 
proach  is  a  relatively  complicated  scanning  mechanism, 
compared  to  seeker  heads  that  scan  in  azimuth. 
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17.0  CONCLUSIONS 


A  signal  processing  algorithm  (MRSI)  has  been 
devised  to  maximize  the  detection  range  of  imaging 
IR  seekers  against  ship  targets.  The  algorithm  com¬ 
prises  a  bank  of  two-dimensional  spatial  filters,  hav¬ 
ing  rectangular  kernels  matched  to  targets  of  various 
aspect  ratios  and  sizes,  implemented  in  a  separable 
form  that  appears  well  suited  to  processing  in  real-time 
image  data  acquired  sequentially  by  column.  The  per¬ 
formance  of  MRSI  has  been  studied  analytically  for 
simple  bar-target-in-noise  images  and  studied 
numerically  via  simulations  performed  on  a  number 
of  measured  and  simulated  ship  images. 

The  SNR  gain  provided  by  MRSI  (relative  to  hot¬ 
spot  detection)  can  be  estimated  as 

where  NP1X  is  the  number  of  pixels  occupied  by  the 
target  in  the  image,  and  =  0.81,  ±10%.  The  pro¬ 
cessor  efficiency,  ij,  is  remarkably  invariant  with  respect 
to  sensor/ship  range,  viewing  aspect,  and  ship  class. 


The  principal  simulation  results  are  plots  of  detec¬ 
tion  probability  (PD)  vs.  SNR  (Figs.  42  and  43),  bas¬ 
ed  on  2300  image-based  Monte  Carlo  trials.  A  curve 
fit  to  the  numerically-developed  probabilities  of  detec¬ 
tion  appears  to  provide  a  universal  curve  that  can  be 
used  predictively  for  ships  of  differing  class,  observ¬ 
ed  at  a  variety  of  aspects  and  distances.  Interestingly, 
the  curve  fit  to  MRSI/Monte  Carlo  PD  results  coin¬ 
cides  with  a  previously  published  curve  fit  to  data  ob¬ 
tained  from  visual  detection  experiments  performed 
with  human  observers  of  TV  displays  [3]. 

Section  15.0  describes  two  ASM  system  concepts, 
distinguished  by  the  means  of  target  classification,  viz., 
classification  performed  on-board  the  missile  by  an 
ATR  algorithm  or  off-board  by  a  person  inspecting 
imagery  telemetered  from  the  missile.  Using  MRSI  to 
cue  a  slow-scan/high  sensitivity  “classification  mode” 
results  in  substantial  predicted  improvements  for  both 
system  concepts. 

In  Section  16.0  we  present  a  new  method  for  sup¬ 
pressing  false  alarms  caused  by  cloud  reflections  from 
the  sea  surface. 
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PROGRAM  MRSIPC 

C  29  APRIL  1987  12:40  PM 

REAL*8  DSEED 

REALM  GAIN 

COMMON  /SEED/  DSEED ,NSEEDS , ISKIP 

COMMON  /XINIT/  KK1 ,KK2 ,K(5),IM(5) 

COMMON  /GAYN/  GAIN (1000) 

DATA  RMS , BCKGND  /l. 0,100.0/ 

C 

CALL  RDINPT 
C 

NTNRS  =  5 

DO  400  ITNR=1, NTNRS 
I TRIAL  =  0 

TNR  =  0.05  +  FLOAT (ITNR)  *  0.10 
TGT  =  BCKGND  +  TNR  *  RMS 
C 

DO  100  ISEED=1 ,NSEEDS 
ITRIAL  =  I TRIAL  +  1 
DSEED  =  DSEED  +  1.0D-03 
C 

C  CREATE  NOISE  FRAME... 

C 

CALL  GAUS (RMS) 

CALL  ADNOIS (0 , 1 . ) 

C 

C  CALIBRATE  MRSI  USING  NOISE-ONLY  FRAME... 

C 

ITST  =  ISKIP  *  (ISEED  /  ISKIP) 

IF(ITST.EQ.ISEED)  WRITE(22 , 210) 

CALL  INIT (KK2P1 ) 

DO  200  KK=1 , KK2P1 
CALL  XAV(KK) 

IMKK  =  IM(KK) 

DO  200  I  =  1 , IMKK 
CALL  IAV(I ,KK) 

CALL  BES(I ,  KK) 

CALL  NE1 (I ,KK) 

200  CONTINUE 

CALL  NE2 

C  CALL  SI GOUT 

C 

C  CALIBRATION  PROCESS  IS  COMPLETE. 

C  ADD  SHIP  IMAGE  TO  NOISE  FRAME... 

C 

CALL  ADSHIP( TGT, BCKGND) 

CALL  ADNOIS (1,1.) 

C 

C  USE  MRSI  TO  LOOK  FOR  SHIP... 

C 

CALL  INIT(KK2P1 ) 

DO  300  KK  =  1.KK2P1 
CALL  XAV(KK) 

IMKK  =  IM(KK) 

DO  300  I  =  1 , IMKK 
CALL  IAV(I ,KK) 

CALL  BES(I , KK) 

CALL  TD(I ,KK) 

300  CONTINUE 

CALL  TDOUT(TNR,GAIN(ITRIAL), ISEED) 


non 
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100  CONTINUE 

IF (ITRIAL.LT. I SKIP)  GO  TO  400 
GENERATE  SUMMARY  OF  SNR  GAIN  PERFORMANCE... 

CALL  GOUT(ITRIAL, I5KIP) 

400  CONTINUE 

210  FORMAT (5X, 'BEGIN  CALIBRATION  PROCESS') 

STOP 

END 

C 

£**************** **** * A****** ********************* ***** A* ****** A******** 

C 

SUBROUTINE  RDINPT 
C 

C**  READS  INPUT  DATA 
C 

LOGICAL*l  TEMP (64, 2 56) 

LOGICALM  LSHIP(64,256) 

INTEGERM  ITEMP4(64 , 256) ,  ISKIP 

REALM  ARAY ,  SHIP 

REAL *8  DSEED 

COMMON  /SEED/  DSEED ,NSEEDS , ISKIP 

COMMON  /XPARAM/  IMAX , KMIN , KMAX 

COMMON  /SPARAM/  JMAX ,NMAX , JAVG , JDEC 

COMMON  /RAWFRM/  ARAY(64 , 256 ) , SHIP (64 , 256 ) 

COMMON  /FILTER/  TS , FCO , NORDER , ISET 
EQUIVALENCE  ( LSHIP , ITEMP4 ) 

C 

READ(21 , 105)  DSEED, NSEEDS, ISKIP 
WRITE (23 , 106)  DSEED , NSEEDS , ISKIP 

105  F0RMAT(1X,D14 .7 , IX, 13, IX. 12) 

106  FORMAT (IX, 1 DSEED= ' ,D14.7,3X, ' NSEEDS= ' ,I3,3X, 'ISKIP=' ,12,/) 

C 

DO  50  I  =  1,64 

50  READ(20 , 100)  (TEMP(I,J)  ,  J=l,256) 

C 

DO  60  J=1 , 256 

DO  60  1=1,64 

LSHIP (I , J)  =  TEMP ( I , J ) 

SHIP (I , J)  =  ITEMP4(I , J)  /  255 
60  CONTINUE 
C 

READ (21,110)  IMAX , JMAX , KMIN , KMAX , JAVG , JDEC , NMAX 
WRITE (23,120)  IMAX , JMAX , KMIN , KMAX , JAVG , JDEC , NMAX 
READ (21,115)  TS, FCO, NORDER, ISET 
WRITE (23 , 125)  TS , FCO , NORDER , ISET 
100  FORMAT ( 2 ( 1 28A1 ) ) 

110  FORMAT(2 (13 ,1X) ,4(12, IX) ,13) 

115  FORMAT(F9 .6 , IX , F7 . 1 , IX , 12 , IX , 12 ) 

120  FORMAT (IX, 'IMAX=' , 13 ,2X, ' JMAX= 1 , 13 , 2X, ' KMIN= 1 , 12 ,2X, ' KMAX= ' ,I2,2X, 
+  ' JAVG=' ,I2,2X, ' JDEC*' ,I2,2X, ' NMAX= ' ,13,/) 

125  FORMAT (IX, 'TS=' ,F9 .6 , 3X, ' FC0= ' ,F7.1,3X, 'NORDER*' ,I2,3X, 'ISET=' , 

+  12,/) 

RETURN 

END 

C 

C**** ******* ***************************************** ******************* 

c 
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C 

C 


INPUT  (FROM  RDINPT): 

ARAY(T, J) 

IMAX , JMAX 
KMIN , KMAX 
JAVG 
JDEC 
NMAX 


1=1,2 . IMAX 

J=1 ,2 . JMAX 

MIN/HAX  VALUES  (CROSS-SCAN  RESOLUTION  PARMS.) 

IN-SCAN  AVERAGING  FACTOR 

IN-SCAN  DECIMATION  FACTOR 

MAX.  VALUE  (IN-SCAN  RESOLUTION  PARAMETER) 


OUTPUT  /  DATA  FRAME  WITH  IN-SCAN  DECIMATION: 

V(I,J)  J=1 , 2 . JMAX2 

OUTPUT  /  CROSS-SCAN  AVERAGING  PARAMETERS: 

KK  INDEX  FOR  K(KK) ,IM(KK) 

=  1,2,3,..., (KK2+1) 

K(KK)  CROSS-SCAN  RESOLUTION  PARAMETER 

=  1,2, 4, 8 . KMAX  (POWERS  OF  2) 

IM(KK)  MAX.  VALUE  OF  I  FOR  V(I,J) 

KK1  SMALLEST  VALUE  OF  KK  FOR  WHICH  CROSS-SCAN 

PROCESSING  SHOULD  BE  PERFORMED 

KK2  NO.  OF  CROSS -SCAN  AVERAGING  STAGES 

OUTPUT  /  IN-SCAN  AVERAGING  PARAMETERS: 


NN 

NN2 

N(NN) 

JMAX2 


INDEX  FOR  N(NN) 

=1,2,3, ... , (NN2+1 ) 

NO.  OF  PARALLEL  CHANNELS  OF  IN-SCAN  PROCESSING 

IN-SCAN  RESOLUTION  PARAMETER 
=1,2, 4, 8 . NMAX  (POWERS  OF  2) 

MAX.  VALUE  OF  J  FOR  V(I,J) 


INPUT  COMMON  BLOCKS  (FROM  RDINPT) 

COMMON  /XPARAM/  IMAX, KMIN, KMAX 
COMMON  /SPARAM/  JMAX, NMAX, JAVG, JDEC 
COMMON  /RAWFRM/  ARAY(64,256) 

OUTPUT  COMMON  BLOCKS 

COMMON  /XINIT/  KK1 ,KK2 ,K(5) , IM(5) 
COMMON  /SINIT/  NN2 ,N(7 ) , JMAX2 
COMMON  /NUFRMS/  V(64 , 256 ) , VIKK(256 , 7 ) 


KK1  =  INT(1 .443*ALOG(FLOAT(KMIN) ) ) 
KK2  =  INT( 1 .443*ALOG(FLOAT(KMAX) ) ) 

KK2P1=KK2+1 
DO  100  KK=1 , KK2P1 
K(KK)  =  2**(KK-1) 


71 


PAGE  004 

LISTING  OF  F1F.SRA.MRSIPC  29  APR  87 


IM(KK)  =  I MAX  -  K(KK)  +  1 
100  CONTINUE 
C 

NN2  =  INT(1 .443*AL0G(FL0AT(NMAX) ) ) 

C 

NN2P1=NN2+1 
DO  200  NN=i,NN2Pl 
N (NN)  =  2** (NN-1 ) 

200  CONTINUE 
C 

JAVGM1  =  JAVG  -  1 
JMAX1  =  JMAX  -  JAVGM1 
JMAX2  =  JMAX1/JDEC 
C 

XJAVG  =  FLOAT (JAVG) 

JDECM1  =  JDEC  -  1 
C 

DO  300  1=1 , IMAX 
C 

C  PERFORM  IN  SCAN  AVERAGING  . . . 

C 

DO  310  J=1 , JMAX1 
SUM  =  0.0 
JPJ=J+ JAVGM1 

DO  320  JSUM=J,J?J 

SUM  =  SUM  +  ARAY ( I , J3UM) 

320  CONTINUE 

V(I,J)  =  SUM/XJAVG 
310  CONTINUE 

C 

C  PERFORM  IN-SCAN  DECIMATION  . . . 

C 

JHOP  =  - JDECM1 
DO  330  J=1 , JMAX2 
JHOP  =  JHOP  +  JDEC 
V(I,J)  =  V( I , JHOP) 

330  CONTINUE 

300  CONTINUE 
RETURN 
END 
C 

C*********************************************************************** 

c 

SUBROUTINE  XAV(KK) 

C 

C**  CROSS-SCAN  AVERAGING  PROCESSOR,  KK-SECTION 
C 

COMMON  /XINIT/  KK1 , KK2 , K ( 5 ) , IM ( 5 ) 

COMMON  /SINIT/  NN2 ,N(7 ) , JMAX2 
COMMON  /NUFRMS/  V(64 , 256 ) , VIKK( 256 , 7 ) 

COMMON  /WORK/  TEMPI (64) ,TEMP2 ( 256 , 7 ) 

C 

IF (KK.EQ . 1 )  GO  TO  100 
IMKK  =  IM(KK) 

KKM1  =  KK  -  1 
DO  200  J=1 , JMAX2 
C 

DO  300  1=1, IMKK 
IPKM1  =  I  +  K(KKM1) 

TEMPI (I)  =  (V(I,J)  +  V(IPKM1 , J) )  /  2.0 
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* 


300  CONTINUE 

C 

DO  400  1=1 , IMKK 
V(I,J)  =  TEMPI (I) 

400  CONTINUE 

200  CONTINUE 
100  CONTINUE 
RETURN 
END 
C 

£* ********* ***************** *** ***** ************************************ 
C 

SUBROUTINE  IAV(I,KK) 

C 

C**  IN-SCAN  AVERAGING  PROCESSOR  /  I-K  SECTION 
C 

COMMON  /SINIT/  NN2 ,N(7 ) , JMAX2 
COMMON  /NUFRMS/  V(64,256) ,VIKK(256 ,7) 

COMMON  /WORK/  TEMPI (64) ,TEMP2(256 ,7) 

C 

DO  50  J  =  1 , JMAX2 
VIKK( J , 1 )  =  V(I,J) 

50  CONTINUE 
C 

C  RECURSIVE  FILTER  INITIALIZATION:  VIKK(N(NN) ,NN) 

C 

WNNM1  =  V(I,1) 

NN2P1=NN2+1 
DO  100  NN=2 ,NN2P1 
SUM  =  0.0 
J1  =  1  +  N(NN-l) 

J2  =  N(NN) 

C 

DO  200  J=J1,J2 
SUM  =  SUM  +  V(I,J) 

200  CONTINUE 

C 

WNN  =  WNNM1  +  SUM 
VIKK(J2,NN)  =  WNN/FLOAT(J2) 

WNNM1  =  WNN 
100  CONTINUE 
C 

C  RECURSIVE  FILTER 
C 

JM2M1=JMAX2-1 
DO  300  NN=2 ,NN2P1 
NNN=N(NN) 

C 

DO  300  J=NNI'I, JM2M1 
JP1  =  J+l 
JP1MN=JP1 -NNN 

SUM  =  V(I,JP1)  -  V(I , JP1MN) 

SUM  =  SUM/ FLOAT (N(NN)) 

VIKK(JPl.NN)  =  VIKK(J.NN)  +  SUM 
300  CONTINUE 
C 

RETURN 

END 

C 

£******************************* **************************************** 
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C 

SUBROUTINE  ADSHIP (TGT , BCKGND ) 

C 

COMMON  /XPARAM/  IMAX , KMIN , KMAX 
COMMON  /SPARAM/  JMAX.NMAX, JAVG, JDEC 
COMMON  /RAWFRM/  ARAY(64 , 256 ) .SHIP (64 , 2S6) 

C 

DELTGT  =  TGT  -  BCKGND 
C 

DO  10  J=1 , JMAX 
DO  10  1=1, IMAX 

ARAY ( I , J )  =  BCKGND  +  DELTGT*SHIP (I , J ) 

10  CONTINUE 
C 

RETURN 

END 

C 

£***********★*****************************★****************************★ 

C 

SUBROUTINE  GAUS(SIG) 

C 

C**  CREATES  A  FRAME  OF  GAUSSIAN  NOISE  WITH  A  WHITE  OR  DIGITAL 
C  BUTTERWORTH  POWER  SPECTRUM. 

C 

C  INPUTS  REQUIRED: 

C  SIG  =  RMS  VALUE  OF  OUTPUT  FRAME 

C  DSEED  =  SEED  FOR  RANDOM  NUMBER  GENERATOR 

C  ISET  =  0  ,  WHITE  NOISE 

C  1  ,  DIGITAL  BUTTERWORTH  FILTER  (STEARNS) 

C 

C  IF  (ISET  .EQ.  1)  ,  ALSO  NEED: 

C 

C  TS  =  CAMPLING  INTERVAL 

C  FCO  =  3  DB  CUT-OFF  FREQUENCY 

C  NORDER  =  FILTER  ORDER 

C 
C 

C  OUTPUT : 

C  GNOISE(64,256)  =  OUTPUT  DATA  FRAME 

C 

C  NOTE : 

C  1)  DATA  ARE  WHITE  IN  THE  FIRST  COORDINATE  AND  BUTTERWORTH- 
C  FILTERED  IN  THE  SECOND  COORDINATE 

C  2)  ISET=1  OPTION  UNAVAILABLE  ON  PC  VERSION 

C 

C  REALM  X(16384)  ,Y(16384) 

REALM  DSEED 

COMMON  /SEED/  DSEED ,NSEEDS , ISKIP 
COMMON  /GSNS/  GNOISE(64, 256) 

COMMON  /FILTER/  TS , FCO , NORDER, ISET 
DATA  NR  /16384/ 

C 

C  GENERATE  WHITE  GAUSSIAN  NOISE. 

C 

CC  CALL  GGNML( DSEED, NR, X) 

CALL  GGNML$( DSEED, NR) 

CC 

CC  IF(ISET  ,NE.  0)  GOTO  50 

CC 

CC  DO  25  K=1 ,NR 


* 

■ 
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CC  Y(K)  =  X(K) 

CC25  CONTINUE 
C  GOTO  60 

C 

C  APPLY  BUTTERWORTH  FILTER. 

C 

C  50  CALL  LOBWF(X, Y,NR,TS , FCO ,NORDER ,5,1) 

C 

C  CALCULATE  STATISTICS  OF  FILTERED  DATA. 

C 

C60  CALL  STAT$$(Y,NR,YMIN,YMAX,YAVE, YVAR.YSIGMA) 

CALL  STATS  $ (GNOI SE , NR , YMIN , YMAX , YAVE , YVAR , YS IGMA , A , B ) 

C 

C  NORMALIZE  THE  OUTPUT  DATA  TO  ZERO  MEAN  AND  STD. DEV.  =  'SIG.' 

C 

C  =  SIG/YSIGMA 
DO  100  J=1 ,256 
DO  100  1=1,64 
C  N  =  J+(I-1)*256 

C  GN0ISE(I , J)  =  C*(  Y(N)  -  YAVE) 

GN0ISE(I , J)  =  C*(GNOISE(I , J)  -  YAVE) 

100  CONTINUE 
RETURN 
END 
C 

£*****★**** ***** AAA *★>•*** ****★★★* ****** * ** ***************** ************* 


c 

c 


c 


10 


SUBROUTINE  ADNOIS(IOPT.TNR) 


COMMON  /XPARAM/ 
COMMON  /SPARAM/ 
COMMON  /RAWFRM/ 
COMMON  /GSNS/ 


I MAX , KMIN , KMAX 
JMAX.NMAX, JAVG.JDEC 
ARAY(64 ,256) 
GNOISE(64 , 256) 


OPT  =  FLOAT (I OPT) 

DO  10  J=1 , JMAX 
DO  10  1=1 , IMAX 

ARAY ( I , J )  =  OPT  *  ARAY ( I , J )  +  GNOISE(I, J)/TNR 

CONTINUE 

RETURN 

END 


C 

£********************  ******  ***********************  *******  *  *******  ******* 


c 


SUBROUTINE  NE1(I,KK) 


NOISE  ESTIMATOR  -  ESTABLISHES  RMS  VALUES  OF  THE  WAVEFORMS  OUTPUT 

BY  1 BES(I ,KK) '  UNDER  CLOSED-COVER  CONDITIONS 


REAL*4 

COMMON  /XINIT/ 
COMMON  /SINIT/ 
COMMON  /NUFRMS/ 
COMMON  /SIGS/ 


X(1536) 

KK1,KK2,K(5),IM(5) 

NN2,N(7) , JMAX2 

V(64 , 256 ) ,VIKK(256 ,7 ) 

5IG( 5 , 7 ) , XSUM(64 ,5,7) ,XSQR(64 ,5,7) ,XNUM(64 ,5,7) 


NN2P1  =  NN2+1 
DO  10  NN=1 ,NN2P1 
INN  =  1/NN 

JSTART  =  2*N(NN)+INN 
JSTOP  =  JMAX2-N(NN)-INN 
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DO  20  J=J START, J STOP 

NI  =  NI+1 

X(NI )  =  VIKK( J ,NN) 

CONTINUE 

CALL  STAT$$(X,NI,A,B,C,D,E,XSUM(I,KK,NN) ,XSQR(I ,KK,NN) ) 
XNUH ( I , KK , NN )  =  FLOAT (NI) 

CONTINUE 


RETURN 

END 

C 

j*********************************************************************** 

C 

SUBROUTINE  NE2 
C 

C**  NOISE  ESTIMATOR  (CONTINUATION  OF  'NEl 1 ) 

C 

COMMON  /XINIT/  KK1 ,KK2 ,K(5) , IM(5) 

COMMON  /SINIT/  NN2 ,N(7) , JMAX2 

COMMON  /SIGS/  SIG(5 ,7 ) ,XSUM(64 ,5,7), XSQR(64 ,5,7) , XNUM(64 ,5,7) 


KK1P1  =  KK1+1 
KK2P1  =  KK2+1 
NN2P1  =  NN2+1 

DO  100  KK  =  KK1P1 , KK2P1 
DO  100  NN  =  1.NN2P1 
XAVE  =0.0 

TSQR  =0.0 

RTOT  =0.0 

SIG(KK,NN)  =  0.0 
IMKK  =  IM(KK) 

DO  200  I  =  1 , IMKK 
XAVE  =  XAVE  +  XSUM(I ,KK ,NN) 
TSQR  =  TSQR  +  XSQR(I ,KK,NN) 
RTOT  =  RTOT  +  XI'IUM  ( I ,  KK ,  NN ) 
200  CONTINUE 

IF(RTOT.LE.O.O)  GO  TO  100 

XAVE  =  XAVE  /  RTOT 

ASQR  =  TSQR  /  RTOT 

VAR  =  ASQR  -  XAVE  *  XAVE 

VAR  =  VAR  *  RTOT  /  (RTOT  -  1.0) 

SIG(KK,NN)  =  SQRT(VAR) 

100  CONTINUE 


C  THIS  IS  A  SHORTCUT,  PROVIDING  QUICK  ANSWERS  FOR  WHITE  NOISE... 
C 

CC  RT2  =  SQRT(2 . ) 

CC  Cl  =  0.0 

CC  DO  300  NN  =  3.NN2P1 

CC  Cl  =  Cl  +  1.0 

CC  C2  =  RT2**C1 

CC  DO  300  KK  =  KK1P1.KK2P1 

CC  SIG(KK.NN)  =  SIG(KK,2)/C2 

CC300  CONTINUE 
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C 

RETURN 

END 

C 

*★★★*******★★  **  A**  *********  A  *  A  A  **********  A**  A**  *********  A  A*  A  A  A'k'k'k'k 

c 

SUBROUTINE  SIGOUT 
C 

C**  PRINTS  OUTPUT  OF  SUBROUTINE  'NE'. 

C 

COMMON  /XINIT/  KK1 , KK2 ,K(5),IM(S) 

COMMON  /SINIT/  NN2 ,N(7 ) , JMAX2 

COMMON  /SIGS/  SIG(5 ,7 ) ,XSUM(64 ,5,7) , XSQR(64 ,5,7) ,XNUM(64 ,5,7) 

C 

KK1P1  =  KK1+1 
KK2P1  =  KK2+1 
NN2P1  =  NN2+1 
WRITE(23 , 105) 

WRITE( 23 ,110)  ( (SIG(KK,NN)  ,  NN=1,NN2P1)  ,  KK=KK1P1 , KK2P1 ) 

105  FORMAT (IX,/, IX. 'SIGMAS' ) 

110  FORMAT (6(1X,F7.3)) 

C 

RETURN 

END 

C 

*********  *********  ***  *****************  *************  ***  A*  ***  *  *************** 

c 

SUBROUTINE  BES(I,KK) 

C 

C**  BACKGROUND  ESTIMATION  &  SUBTRACTION 
C 

COMMON  /XINIT/  KK1 , KK2 ,K^5),IM(5) 

COMMON  /SINIT/  NN2 . N ( 7 ) , JMAX2 
COMMON  /NUFRMS/  V(64 , 256 ) , VIKK( 256 , 7 ) 

COMMON  /SIGS/  SIG( 5 , 7 ) , X5UM(64 ,5,7) ,XSQR(64, 5 , 7 ) ,XNUM(64 ,5,7) 

COMMON  /WORK/  TEMPI (64) ,TEMP2(256 ,7) 

C 

NN2P1  =  NN2+1 
C 

DO  200  NN=1 , NN2P1 
INN  =  1/NN 

JSTART  =  2*N(NN)+INN 
JSTOP  =  JMAX2-N (NN ) - INN 
NNM1  =  NN-l+INN 
C 

DO  200  J= JSTART, JSTOP 
JL  =  J-v3*N(NN))/2  -  INN 
JR  =  J+N(NN)+INN 

VHAT  =  (VIKK(JL.NNMl)  +  VIKK( JR ,NNM1 ) )/2.0 
TEMP2(J ,NN)  =  VIKK(J.NN)  -  VHAT 
200  CONTINUE 
C 

DO  300  NN=1 , NN2P1 
INN  =  1/NN 

JSTART  =  2*N(NN)+INN 
JSTOP  =  JMAX2-N(NN)-INN 
C 

DO  300  J» JSTART , JSTOP 
VIKK(J.NN)  =  TEMP2( J ,NN) 

300  CONTINUE 
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C 

RETURN 

END 

C 

*  ***  ***:★★**  ******  *****  A  *  ****  k'k'k  A  *********  *****  ***********  ********  A********* 

c 

SUBROUTINE  TD(I ,KK) 

C 

C**  THRESHOLD  /  DECISION  PROCESSOR 
C 

COMMON  /XINIT/  KK1 , KK2 , K( 5 ) . IM( 5) 

COMMON  /SINIT/  NN2 ,N(7 ) , JMAX2 
COMMON  /NUFRMS/  V(64 , 256) ,VIKK(256 , 7 ) 

COMMON  /SIGS/  SIG( 5 ,7 ) , XSUM( 64 ,5,7), XSQR(64 ,5,7) , XNUM(64 ,5,7) 

COMMON  /XCDNS/  JX(64 . 5 ) ,NNX(64 , 5) , VX(64 , 5) 

C 

NN2P1  =  NN2+1 
C 

DO  100  NN=1,NN2P1 
INN  =  1/NN 

JSTART  =  2*N(NN)+INN 
JSTOP  =  JMAX2-N(NN)-INN 
C 

DO  100  J= JSTART, JSTOP 

VIKK( J ,NN)  =  VIKK(J,NN)/SIG(KK,NN) 

100  CONTINUE 
C 

JX(I,KK)  =  1 

NNX (I , KK )  =  1 

VX ( I , KK )  =  VIKK(3 , 1 ) 

C 

DO  400  NN=1,NN2P1 
INN  =  1/NN 

JSTART  =  2*N(NN)+INN 
JSTOP  =  JMAX2-N(NN) -INN 
C 

DO  400  J= JSTART, JSTOP 

IF(VX(I ,KK) .GT.VIKK( J ,NN) )  GOTO  400 

JX ( I , KK )  =  J 

NNX ( I , KK )  =  NN 

VX ( I , KK )  =  VIKK( J ,NN) 

400  CONTINUE 
C 

RETURN 

END 

C 

£********* A******************* **************** A*** ********************** 

C 

SUBROUTINE  TDOUT(TNR,GAIN, ISEED) 

C 

C**  PRINTS  OUTPUT  OF  SUBROUTINE  ' TD * 

C 

COMMON  /XCDNS/  JX(64, 5) ,NNX(64 , 5) , VX(64 , 5) 

COMMON  /SPARAM/  JMAX , NMAX , JAVG , JDEC 
C 

IF(ISEED.EQ.l)  WRITE(23 ,95) 

IF(ISEED.EQ.l)  WRITE(23 ,100) 

IF(ISEED.EQ.l)  WRITE(23 ,110) 

95  FORMAT (IX, '  ',/) 

100  FORMAT ( '  '.IX,'  TNR  IX  JX  K  N  GAIN') 
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110  FORMAT (2X,  1  .  .  . ') 

210  FORMAT (4X, F5 . 2 , 6X , 12 , 4X , 13 , 4X , 12 , 6X , 12 , 5X , F7 . 3 , / ) 

C 

IX  =  1 

JXO  =  JX(1,1) 

KKX  =  1 

NNXO  =  NNX (1,1) 

VXO  *  VX(1,1) 

DO  300  I  =  1,64 
DO  300  KK  =  1,5 
IF(VX(I,KK).LE.VXO)  GO  TO  300 
IX  =  I 

JXO  =  JX(I ,KK) 

KKX  =  KK 

NNXO  =  NNX (I ,KK) 

VXO  =  VX(I ,KK) 

300  CONTINUE 
C 

J  =  JDEC*JX0 
K  =  2**(KKX-1) 

N  =  JDEC*2**(NNXO-l ) 

GAIN  =  VXO/TNR 

WRITE(23 ,210)  TNR , IX  . J ,K,N .GAIN 
C 

RETURN 

END 

C 

************************************* ****************** *********** 

c 

SUBROUTINE  GOUT ( ITRIAL , I SKIP ) 

C 

C**  PRINTS  OUTPUT  SUMMARY  OF  SNR  GAIN  PERFORMANCE 
C 

RFAL*4  GAIN 

COMMON  /GAYN/  GAIN(IOOO) 

C 

WRITE(23 , 100) 

WRITE (23, 200) 

WRITE(23 ,300) 

WRITE (23, 400) 

WRITE (23, 500) 

100  FORMAT (IX, '  ',//, 

i  *******************************************************  t//) 

200  FORMAT (12X, 'SUMMARY  OF  SNR  GAIN  PERFORMANCE' ) 

300  FORMAT  (12X,  ' . ',//) 

400  FORMAT ( '  ' ,4X, 'TRIALS' ,4X, 'GMIN' ,7X, 'GMAX' ,7X, 'GAVE' , 7X, 'GSIGMA* ) 

500  FORMAT  ( '  '  ,4X,  ' . '  ,3X,  ' . '  ,3(4X,  ' . ')) 

600  FORMAT (6X, I3,4(4X,F7.3)) 

C 

DO  700  ITRL  =  2, ITRIAL 
ITST  =  ISKIP  *  (ITRL  /  ISKIP) 

IF(ITST.NE.ITRL)  GO  TO  700 

CALL  STATS  $ ( GAIN , ITRL , GMIN , GMAX , GAVE , GVAR , GSIGMA, A, B) 

WRITE (23,600)  ITRL . GMIN , GMAX , GAVE .GSIGMA 
700  CONTINUE 
C 

WRITE(23 , 100) 

RETURN 

END 

C 
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DO  50  IVEC*1 ,NR 

GENERATE  NEW  PAIR 
GENERATE  FIRST  UNIFORM  NUMBER 
IF  (K.EQ.2)  GO  TO  40 

J  *  IFIX(SNGL(Z)) 

C  X  =  Z  -  DFLOAT(J)  +  DFLOAT(J)/P 

X  =  Z  -  DBLE(J)  +  DBLE ( J)/P 

XI  =  TWO  *  X  -  ONE 
IF(Xl.EQ.O)  GO  TO  10 
IF(Xl.GT.ONE)  GO  TO  10 

C  GENERATE  SECOND  UNIFORM  NUMBER 

20  Z  *  IA  *  AX 

J  *  IFIX(SNGL(Z) ) 

C  AX  =  Z  -  DFLOAT(J)  +  DFLOAT(J)/P 

AX  =  Z  -  DBLE(J)  +  DBLE ( J)/P 

X2  =  TWO  *  AX  -ONE 
IF(X2.EQ.O)  GO  TO  20 
IF(X2.GT.ONE)  GO  TO  20 
C  DETERMINE  ACCEPTANCE/REJECTION 

Y1  =  X1*X1  +  X2*X2 
IF(Yl.EQ.O)  GO  TO  10 
IF(Yl.LT.ONE)  GO  TO  30 
GO  TO  10 

C  TRANSFORMATION  TO  NORMAL  DISTRIBUTION 

30  Y1  =  XI  *  DSQRT ( -TWO  *  DL0G(Y1)/Y1) 

Y2  =  Y1  *  (X2/X1) 

XVEC(IVEC)  *  Y1 
GO  TO  50 
40  K  =  1 

XVEC(IVEC)  =  Y2 
C 

50  CONTINUE 
RETURN 
END 
C 

BLOCK  DATA 

INTEGER*2  K 

REAL*8  I , IA , P ,Y2 

COMMON/HOLD/ I , I A , P , Y2 , K 

DATA  K/l/, 1/4194305. DO/, IA/2097153. DO/ 

DATA  P/2147483647. DO/ ,Y2/0. DO/ 

END 
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